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ABSTRACT 


A microcomputer-based, low-cost Omega sensor processor (OSP) has been 
designed, implemented with minimum hardware, and tested. The feasibility of 
implementing a minimum hardware/maximum software system of less complexity 
and lower cost than is currently available has been verified. The design approach 
and hardware/software details are presented for this low-cost navigation aid. 

The long-range project goal of developing a general aviation Omega 
navigation system (ONS) for a cost under SlOOO is now feasible. The micro- 
computer-based OSP presented here represents neatly all the system hardware 
needed for an airborne ONS. An area navigation routine is the additional 
component needed to complete the ONS for general aviation. This navigation 
routine can readily be added to the present system since it is additional soft- 
ware and not hardware. 

This paper describes the OSP design and engineering approach used In 
achieving this design, m reflecting a structured program design approach, the 
project goal of designing a low-cost OSP is presented initially. Succeeding 
chapters document the efforts made in support of this goal. The OSP operations 
are functionally described in the text, as Is the system hardware and software. 
Circuit details of the hardware modules and program listings for the FORTRAN 
and microcomputer software are reserved for the appendices. 
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INTRODUCTION 


This project was sponsored by the NASA/Langley Research Center under a 
grant effort known as the NASA Tri- University Program in Air Transportation 
Systems. The program was begun in 1971 and includes MIT/ Ohio University/ 
and Princeton University in a cooperative effort to improve the safety/ reli- 
ability/ and efficiency of our National Airspace System. Eacn University offers 
a unique area of specialized interest and expertise, thereby complementing each 
others* efforts; Ohio University's area of specialization being avionics. 

Early in the program's development it was determined that the Tri- 
University Program might make a contribution to general aviation by developing 
a low-cost area navigation system for small aircraft. Studies pointed to the 
Omega navigation system as a low-cost system which could provide all-altitude 
signal coverage, and the major emphasis in the program became the development 
of a lov-cost Omega navigation system for the general aviation pilot. 

While MIT studied the systems integration aspects and Princeton Univer- 
sity the flight dynamics aspects, Ohio University turned its efforts towards 
designing a more sensitive and lower cost Omega receiver than was commercially 
available. An early result was an Omega sensor processor incorporating all 
TTL circuitry, and this unit was flight- tested aboard the Avionics Engineering 
Center's DC-3 flying laboratory enroute to Langley Field, Virginia in October, 
1974, After this successful flight demonstration, Ohio University embarked on 
a project to build five Omega prototype receivers for the other Universities 
and NASA. These prototypes were designed using low power CMOS circuitry 
and were delivered to the recipients in September, 1976. 


The Omega sensor processor outputs. Omega phase and phase differences, 
are the inputs needed hy a navigation processor to derive pilot-usable information 
such as cross-track error and down-track distance. Due to the recent availability 
of iow-cost microcomputers, efforts were channeled toward using a commercial ly- 
avai table microcomputer as the navigation processor for the airborne Omega 
navigation system (ONS). When incorporating a microcomputer into the ONS, 
however, it became obvious that, in addition to performing the area navigation 
functions, the sensor processor functions could also be performed in micro- 
computer software f 

Therefore, a project was begun to explore the feasibility of a minimum 
hardware/maximum software OSP. This paper reports results of this project: 
a minimum hardware system has been designed, and the OSP functions of auto- 
matic synchronization, phase tracking, and line-of-position (LOP) generation 
have been implemented in microcomputer software. The microcomputer-based 
OSP thus consists of a 10.2 KHz TRF receiver module, an interface module 
which samples the receiver outputs, and a microcomputer module to process 
the Omega phase data. 

The design approach was to first develop the system hardware (receiver 
and interface modules) and use these modules to supply Omega data to a digital 
magnetic tape recorder. A data base of real Omega ground and flight data was 
collected for use in a FORTRAN language simulation analysis of OSP functions. 
Algorithms to perform the OSP functions were developed and optimized via this 
simulation. The routines were then assembled into microcomputer code and loaded 
into the OSP microcomputer module. With the receiver and interface supplying 
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Omega data to the microcomputer module, the output data consistently equaled or 
exceeded the quality of the all-hardware OSP outputs, thus verifying the feasi- 
bility of the software-based system. 

The hardware and software for the microcomputer-based low-cost OSP 
are functionally described herein, and the type of outputs available and the 
results obtained are presented. The hardware circuit details and software 
routines program listings are documented in the appendices. 

II, OMEGA SENSOR PROCESSOR 

The worldwide Omega navigation system is a network of eight VLF 
transmitting stations. The stations transmit in a time-multiplexed fashion on 
the frequencies of 10.2 KHz, 13.6 KHz, and 11.333 KHz. A brief system 
description is given in Appendix A, while a more thorough treatment is afforded 
by U.S. Coast Guard publication "Omega Navigation System User Handbook". 

C23 

In-depth treatment of Omega system details has been given by Pierce^ and 
Morris. 

Information on the design of airborne receiving and navigating systems 

using the Omega signals is less accessible. Although commercial manufacturers 

are understandably secretive with their receiver hardware/software designs, 

reports on those systems whose development has been funded by government 

agencies should be available to the public (e.g., " Omega Navigation Set: 

An/ARN-99 Final Engineering Report" An all-hardware Omega sensor 

processor designed by Ohio University for NASA is described in "Ohio University 

[53 

Omega Prototype Receiver". 
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An Omega navigation system for airborne use can consist of three func- 
tional units; an Omega sensor processor (OSP), navigation processor, and a cock- 
pit display. The OSP receives the Omega signals and provides amplification and 
signal-to-noise ratio improvement. The OSP must provide Omega data usable 
by the navigation processor. This navigation processor must convert the Omega 
information into pilot-usable information to be displayed on the cockpit display. 

A simple display sufficient for general aviation could take the form of analog 
or digital readout of cross-track error and down-track distance. 

A. Design Approach. The project effort to develop a lov/-cost Omega 
navigation system for genera! aviation has been greatly aided by the recent 
commercial availability of inexpensive microcomputers. The reasons for chosing 
to use a microcomputer to perform navigation processing are well known: high 
reliability, small size, low cost, low power consumption, and off-the-shelf 
availability in addition to software design flexibility. Using a microcomputer 
for sensor orocessing as well as navigation processing has the advantages noted 
above, but it can present some complex systems problems to the designer 
(particularly time sharing of the computer for several functions and t^rioi itizing 
of work). Therefore, this project was undertaken to explore the feasibility of 
implementing an Omega sensor processor with a minimum hardware/maximum 
software combination. 

The design approach taken was to use the previously developed Ohio 
University low-cost front-end receiver module in conjunction with an inter- 
face module to collect a data baseofreal Omega ground and flight data. This 
data base was then used as the input for a FORTRAN language simulation analysis 
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of tl>e OSP funcHons. After the OSP functions had been simulated and optimized 
in FORTRAN, they were rewritten in microcomputer language and tested in the 
final system configurationp The power of tnis approach lies in the fact that the 
software routines were developed not only with computer-generated data as input, 
but they were optimized to perform with the actual Omega data as collected with 
the OSP hardware. 

The OSP is functionally illustrated in Figure 1, and the hardware modules 
and software algorithms are described in the following two sections. For the 
purpose of data collection, the microcomputer module of Figure 1 was replaced 
with a digital magnetic tape recorder. This setup was used tn the laboratory for 
the collection of ground-based data, while the collection of flight data was 
accomplished with this configuration onboard Ohio University’s DC-3 aircraft 
N4002. Complete documentation of the Omega phase data base is given in 
Appendix B. 

Figure 1 shows that the basic system operation is as follows; the receiver 
front-end module receives the Omega signal from the antenna-preamplifier and 
provides filtering, amplification, and hard-limiting before sending Omega zero- 
crossings to the interface module. The microcomputer interface module measures 
and samples the Omega phase information at some predetermined rate and generates 
a microcomputer interrupt signal. The microcomputer (or data recorder) receives 
the digitized phase sample, and this information is used by the OSP software 
routines. 

The OSP software routines for automatic synchronization, phase tracking, 
and LOP generation perform the bulk of the Omega phase processing work. 
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B. System Hardware. The Omega sensor processor design is intended to 
be a minimum hardware configuration with emphasis on low cost. The use of a 
general purpose microcomputer to perform as many OSP functions as possible 
reflects this design objeciive. The OSP is functionally illustrated in the block 
diagram of Figure 1, and a description for each of the hardware modules follows. 
Circuit details for each module are given in Appendix C. 

1 . Receiver. The antenna used for airborne operation is the air- 
craft’s ADF sense antenna. If the aircraft is not equipped with an ADF sense 
wire, another high-impedance E-field antenna may be used (such as a bent dipole 
nav-com antenna). For ground-based operation it is desirable to mount the antenna 
as high above buildings, etc., as possible and in an area free of 60 Hz inter- 
ference {10.2 KHz being the 170th harmonic of 60 Hz). 

The signal developed across the antenna is amplified by a preamplifier 
which is colocated with the antenna. This colocation is necessary to overcome 
signal attenuation due to capacitance of the receiver module lead-in cable. 

The preamplifier is primarily an impedance-matching device and has two outputs. 
The Omega output provides a gain of 20 decibels at 10.2 KHz in a -3 dB band- 
width of several kilohertz. The phase shift is adjusted to zero degrees at 
10.2 KHz. The ADF output provides a flat response of -6 dB across the band 
from 10 KHz to 1 MHz, and the phase shift is zero degrees across the ADF band. 
The preamplifier response characteristics are shov-^n In Figure 2. The preamplifier 
has been described by Wright. 

Accepting the signal from the preamplifier is a 10.2 KHz TRF receiver 
front-end module (see Figure 3). The front-end module provides additional gain. 
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Figure U Summary Block Diagram - Microcomputer-Based Omega Se..sor Processor 
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Figure 2. Preamplifier Characteristics. 
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narrowband filtering, and hard-limiting (i.e. provides pulse outputs coincident 
with the 10.2 KHz Omega zero-crossings). The unit consists of three integrated 
circuit chips and two ceramic filters providing a bandwidth of 30 Hz at 10.2 KHz. 

The module is equipped with an RF gain control which permits gain adjustment 
from 160 to 200 dB. The front-end frequency response curve is shown in Figure 
4. As can be seen from this curve, the response is down nearly 80 dB at the other 
Omega frequencies of 11,333 KHz and 13.6 KHz. The sensitivity and recommended 
operating range is 1 microvolt to 10,000 microvolts. The front-end module has 

been described by Burhans, 

2. Computer Interface . After the bandpass filtering and hard- limiting 
provided by the front-end, pulses coincident with the Omega zero-crossings are 
sent to the microcomputer interface module. This module is also illustrated in 
Figure 3. The functions of phase detection and data sampling are performed by 
this interface hardware, and it also generates the interrupt signals necessary for 
reading the phase data into the microcomputer. Besides the Omega zero-crossings, 
a clock signal of 1.3056 MHz is used by the microcomputer interface module. 

This clock signal is synthesized from a 5 MHz temperature compensated crystal 
oscillator. The choice of 5 MHz as the basic oscillator frequency permits operation 
by an external, more stable clock of this standard frequency should this be desired. 

7 

Phase detection is performed in an open loop fashion: the 2 x 10.2 KHz 
clock signal is input to a 7-bit reference counter, and the seven counter outputs 
are used as data inputs to a 7-bit latch. The Omega zero-crossings clock the 
latch, thus providing a relative phase measurement of Omega compared to the 
local clock approximately 10,200 times per second. 
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Figure 4. Frc 't-’End Module Response Curve. 








There are af least two considerations which suggest sampling the Omega 
phase measurements at a rate slower than the 10,2 KHz carrier rate: 

(1) Since an extremely precise measurement is not required, it 
is unnecessary to use the phase measurement of every cycle of 
the carrier to gain a filtered estimate of the true phase. 

(2) Using the data at the maximum available raw-data rate 

implies a prohibitively high interrupt rate for most micro- 
computers. ' 

The choice of a particular sampling rate is also based on these consider- 
ations. in this implementation a sampling rate of 100 Hz has been chosen for the 
following reasons: 

(1) Sampling of the 15 Hz post-detection bandwidth phase at 100 Hz 
is a rate high enough to avoid severe aliasing of the sampled data 
spectrum, and it provides a sufficient statistical base of data 
points for each Omega time-slot. 

(2) 100 Hz is an acceptable microcomputer interrupt rate, and it 
is easily generated from the existing timing c.icuitry. 

The microcomputer interface module derives the 100 Hz sampling frequency, 
samples the phase data at this rate, and sends an interrupt signal to the micro- 
computer. Having received the interrupt signal, the microcomputer services 
the interrupt request by loading the new phase data from the interface module. 

The inexpensive technique known as memory-mapping is used to accomplish this 
loading of data. The microcomputer address bus is decoded, thereby enabling 
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the data onto the data bus and into the microcomputer (as indicated in Figure 3). 
The interface module and memory-mapping have been previously addressed by 
this author in references 8,9, and 10. 

3. Microcomputer. The microcomputer module consists of a commer- 
cially available microcomputer system made up of a central proce.ssing unit (CPU), 
random access memory (RAM), peripheral interface adapter (PIA), teletypewriter 
interface circuitry, and teletypewriter monitor software in permanent read-only 
storage (ROM) as illustrated in Figure 5. The CPU utilizes a sixteen- line address 
bus, enabling addressing up to 65 kilobytes of memory. An eight-bit bidirectional 
data bus allows for data transfer among the CPU, memory, and input/output 
(I/O) devices. 

The computer's clock is a 1 MHz adjustable RC oscillator or crystal 
oscillator, resulting in a machine cycle time of one microsecond. Most Instructions 
require 2 or 3 machine cycles to execute. CPU registers include one accumulator, 
two index (scratch-pad) registers, a processor status register, a stack pointer, 
and a program counter. All of these are eight-bit registers except the program 
counter which is sixteen bits. Appendix C gives the microcomputer specifics, 
and its operation is explained in detail in the manufacturer's manuaIsP ^ 

Figure 6 shows the simple, inexpensive method for output of the micro- 
computer data, which is the same technique used for data input mentioned 
previously. An output latch is memory-mapped to a particular address by 
decoding the microcomputer address bus lines. The address decoder output 
signal is gated with the microcomputer's "write" signal. This gate's output 
enables a latch which stores the logic levels appearing on the data bus at that 
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Figure 6. Memory-Mapped Data Output Port. 
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time. This data is then fed to a digita!-to-analog (D/A) converter for display 
on a chart recorder or meter. 

C. System Software , 

1. Interrupt Servicing . When the microcomputer interface module 
sends an interrupt signal to the microcomputer, the interrupt service routine 
loads the new data into the computer. Then the software process as flow-charted 
in Figure 7 begins. The basic interrupt service routine serves as a system time- 
keeper and data "traffic cop": it loads in new data as it becomes available and 
hands it off to the proper data-processing subroutine. If the OSP is not in sync 
with the Omega transmission format, synchronization must be achieved first. 
Otherwise, the data is sent to the tracking loop subroutine if it is time-slot 
data, whereas the data is simply ignored if it is from an Omega transmission 
gap, A data output subroutine is executed if it is the correct time for it. 

The automatic synchronization, tracking loop, and data output sub- 
routines called by this interrupt service routine are functionally described in 
the subsequent text material, while the FORTRAN and microcomputer language 
programs are listed and described in Appendices E and F. The interrupt service 
routine requires 150 bytes of memory and performs its time and data management 
in a few tuns of microseconds following each Interrupt request signal. 

2. Automatic Synchronization . The Omega system consists of eight 
stations transmitting on three VLF channels in a fixed, time-multiplexed format. 
As described in Appendix A, a unique pattern is formed by the scheduled length 
(in time) of each station“s transmission which repeats every ten seconds. Thus, 
to use the phase measurements supplied by the microcomputer interface module. 
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Figure 7. Interrupt Service Routine Flow Chart. 
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some means musf be provided for synchronizing the OSP timing with the Omega 
transmission pattern. Only after the OSP is "in sync" can the phase of each station's 
signal be followed by tracking filters to give position-fixing information. 

It should be pointed out that precision sync (better than plus or minus 100 
msec) is unnecessary. A .2 second gap between each transmission burst allows 
for propagation delays from transmitter to receiver and for signal build-up on the 
transmitting antenna and through the narrowband receiver module without signal 
overlap into another station's time-slot. Therefore, sync within plus or minus 100 
msec insures that only one station will be received during each of the eight time- 
slots. 

A description of the automatic synchronization procedure simulated in FORTRAN 

and implemented in the OSP microcomputer follows, and a logic flow diagram appears 

in Figure 8, Although the routine is described here in terms of a single-frequency 

OSP, the technique is applicable to multiple-frequency systems as well. Per- 

« 

forming the routine on multiple channels simultaneously would give the OSP a 
redundancy check. By averaging the results obtained for all channels a best 
estimate of the true sync point could be obtained. 

To establish sync a pattern of station transmission lengths is derived from 
the received data, and the entire frame of data is shifted until the received pattern 
coincides with the unique transmission pattern stored in the computer. The received 
pattern of transmission lengths can be established by determining the amount of 
coherence between consecutive iO msec phase samples. A high level of coherence 
indicates the presence of a strong signal, while the lack of phase coherence (i.e. 
phase jitter) is indicative of a weak signal or noise. The absolute value of phase 
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measurements is meaningless during the sync process; only the phase coherence 
from sample to sample is useful in determining the received transmission pattern 
(i.e. in determining whether each sample is signal or noise). 

An Omega phase measurement is supplied by the microcomputer interface 
module every 10 msec (100 Hz sampling rate). Starting at a random point in 
time (power on time) and taking data for ten seconds (one complete Omega frame) 
yields 1000 phase samples. This is all the information needed to accomplish 
sync. The Omega signal is contaminated with noise, r:ausing a series of measure~ 
ments to result In a dispersion of phase points rather than a constant value. The 
units of centicycles (cec) are commonly used in Omega ;argon. One cec »s 
1/100 of an Omega carrier cycle; at 10.2 KHz a cycle is about 16 nautical 
miles long, making one cec equal to about .16 miles. If the phase of the received 
signal with respect to the OSP clock is near the 0 or 99 cec measurement extremes 
(near the edge of a cycle), the measured value may "bobble" between these 
extremes. To remove this bobble each data point above 49 cec Is "mirrored" 
about the 50% full-scale value (mirrored about 49 cec). Points at x " 99 cec 
are thus reflected to 99-x = 0 cec. Since only the phase cohere, ce is of interest 
and not the absolute value of phase, the mirroring process merely puts the data 
into a better form from which coherence can be determined. 

A simple routine yielding a measure of the point-to-point coherence Is 
to take the difference of each two successive phase points. That is, subtract- 
ing each measurement from the previous yields 1000 difference values for the 
10 second data frame. Each difference is inversely proportional to the coherence 
between the two phase samples. 
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By comparing each difference value to a predetermined "threshold constant", 
a decision can be made as to whether signal or noise exists in each 10 msec sample 
slot. However, it is desirable to first smooth or average the difference values 
over several sample slots before comparison to the threshold constant is made. 

If the smoothed difference is less than the threshold, a "1 " is stored for the sample 
slot; whereas, if the smoothed difference is greater than the threshold, a "0" is 

stored for the slot (low difference high coherence signal present). 

After all 1000 comparisons have been performed, a pattern of 1000 I's and 0*s 
represents the signal and noise sample slots, respectively. The "high-low" or 
"on-off" pattern created by the 1000 binits resembles the Omega transmission 
format shown in Appendix A, the amount of resemblance being proportional to 
each station’s signal-to-noise ratio (SNR). 

Note that the 1000 bits can be stored in 125 eight-bit bytes of memory 
(an amount of storage not prohibitive for simple microcomputer systems). Also 
bear in mind that the entire microcomputer will be dedicated to the auto-sync 
routine during the sync process, since no Omega navigation information is avail- 
able until sync is completed. 

The 1000 sample bits are next compared to the stored Omega pattern in 
a bilevel correlation process. Use is made of the time-slot and transmission gap 
information; whenever the received slot binit is equal to the stored slot binit, 
a correlation counter is incremented by one. After all 1000 slots have been 
compared, the stored pattern is shifted by one slot and the correlation of points 
sequence Is repeated. This shift and correlate iteration is continued until all 
i 000 possible shifted patterns hcve been tested. The number of shifts necessary 
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to obtain the highest bllevei correlation value is used to obtain {!ie starting point 
for the "A" time-^lot (assuming the microcomputer's interrupt service routine 
has been keeping time (counting 1000 Hz interrupts) in a modulo 10-second 
foshion since the auto-sync process began). 

Having found the sync-point in the simulation data. Figure 9 shows the 
original 1000 byte record of phase data plotted ten points per line. The higher 
the SNR of the station being received in each time-slot, the less dispersion 
there is in the phase points plotted. 

The OSP operator's a priori knowledge of which stations are on the air 
and which ones will be well received in his particular geographic area can be 
incorporated into the creation of the unique transmission pattern stored in the 
computer. This facilitates the bilevel correlation process, resulting in a more 
accurate sync point. 

Although it was not used in this version of the OSP, for an additional 
confidence measure the entire sync process can be performed several times and a 
best estimate of the sync point computed. If multiple frequency channels are 
available, the process can be performed on each channel and the sync points 
averaged (or a sync point can be considered acceptable only if it results on two 
of three channels, etc.). Another possibility is to save the second highest 
correlation value as well as the highest, and consider the sync point acceptable 
only if the highest value is better than the second highest by some predetermined 
amount. 

in summary, th« technique presented requires very simple microcomputer 
instructions (compare, subtract, shift), uses very little scratch-pad RAM storage 
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Figure 9, Computer Plot of Raw Omega Data, This is the 30 Hz 
bandwidth phase plotted 10 points per line. 


■23 



(125 bytes), and usually accomplishes sync in about 100 seconds (depending on 
the redundancy constraints imposed). The technique has demonstrated acceptable 
accuracy (within ± 100 msec of true sync) and the microcomputer program Is about 
300 bytes long. 

3. PLL Phase Tracking . After the OSP is in sync with the Omega 

transmissions, the phase of each station's signal can be tracked using correlation 

techniques, and SNR improvement can be made. The correlation detection 

structure known as the memory-aided phase locked loop (MAPLL) has been found 

[ 13 ]. 

to be a successful tracking filter for the Omega signals. Since previous 
experience with hardware versions of the MAPLL had revealed its usefulness 
for Omega tracking, it was desirable to duplicate this configuration in software. 

The first-order version of the software memory— aided phase locked loop 
(SMAPLLl) is functionally illustrated in Figure 10. The loop operation can be 
generally explained as followst the loop's bilevel quantizing phase detector com— 

I 

pares the incoming phase measurement with a filtered loop control word (LCW) 
at the 100 Hz sampling rate. The output of the phase detector is on increment 
command if the incoming phase is ahead of the LCW; it is a decrement command 
if the incoming phase lags behind the LCW. The increment/decrement commands 
are f.ltered in a bidirectional loop counter. Loop filtering is achieved by requir- 
ing the loop counter to receive several correlated signals (counts in the same 
direction) before it affects the LCW. When the loop counter has received 
sufficient correlated up/down commands, it outputs an increment/decrement 
command to the LCW register. This new LCW is in turn compared to the new in- 
coming phase measurements in the phase detector. 
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Figure 10. Firsl--Order Sofhvare Memory-Aided Phase Locked Loop (SMAPLLl) Block Diagram. 







The memory-aiding feature allows a single phase locked loop (PLL) routine 
to be utilized to track all eight Omega stations' phase. At the end of a current 
{•ime-slot the LON is stored in a memory location, and the LCW for the next time- 
slot is retrieved from memory and used to track the new time-slot's phase. 

Atmospheric noise in the VLF-Omega frequency range is known to be caused 
primarily by lightning. ^ The noise is made up of both gaussian and impulsive 
noise. The bilevel quantizing phase detector is thought to discriminate against 
the impulsive noise better than a linear detector (in addition to being simpler 
to implement). Omega signal levels are such that the PLL must be able to track 
signals at a -20 dB SNR in the 30 Hz predetection bandwidth. The loop filter 
(bidirectional counter) provides a narrow post-detection bandwidth of a fraction 
of a Hertz to enable tracking of these low-level signals. However, the increment/ 
decrement type of phase detector and narrowband loop filter severely limit the 
tracking speed of the loop. Increased loop filtering gives SNR improvement but 
decreases the maximum tracking speed of the loop, and vice-versa. 

The PLL must be capable of tracking the motion of a general aviation air- 
craft (nominally 200 knoh velocity). The frequency offset of the receiver's 
local oscillator (TCXO) appears as a velocity offset in phase, thereby adding a 
few more tens of knots to be tracked by the PLL. This speed requirement coupled 
with the narrowband filtering restraint dictates the need for a higher order tracking 
loop (i.e. one that can track a constant velocity offset without sacrificing the 
needed integration). Therefore, the second-order memory-aided PLL (SMAPLL2) 
structure illustrated in Figure 11 has been implemented. 
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Figure 11. Second-Order Software Memory-Aided Phase Locked Loop (SK'APLL2) Block Diagram. 














The operaHon of this PLL is similar to that of the first order loop, but a 
second-order loop counter and a velocity register have been added. When 
the first-order loop counter increments or decrements the LCW, it also increments 
or decrements the second-order loop counter. At the end of the time-slot the 
contents of the second-order loop counter are added to the velocity register. 

The velocity register contents are then scaled by a second-order gain factor 
and used to rate-aid the LCW before its use in the next time-slot measurement 
10 seconds later. It should be noted that the SMAPLL2 exhibits first-order loop 
characteristics during each time-slot, with second-order characteristics (rate- 
aiding) only between time-slots. 

Additionally, a SNR measurement is derived from the phase detector 
output. As shown In Figure 11 , the phase detector compares the incoming phase 
mrasurement with the LCW. This is accomplished by differencing the two numbers. 
This difference is a measure of how "tightly" the loop is locked, which is in turn 
a function of the received SNR, The absolute (unsigned) values of these differences 
are accumulated over all measurements in a time-slot. For a high SNR case the 
accumulation will be small (loop tightly locked), while for a low SNR case the 
sum will be large. 

The SAAAPLL2 parameters of first order bandwidth and second order gain 
have been optimized via FORTRAN simulation as documented in Appendix E. 
Worst-case conditions of 200 knot aircraft velocity and 100 knot TCXO clock 
drift in a radial direction to/from Omega stations, as well as a SNR range of 
-20 dB to +20 dB in the 30 Hz predetection bandwidth have been taken into 
account in the optimization procedure. A first order filter of 16 counts (about 
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.05 Hz bandwidth) combined with a second order gain of .25 have been chosen 
to yield the best combination of transient and steady-state responses. 

Figure 12 shows SMAPLL2 transient (loch-up) response to 180 degree step 
and 300 knot ramp inputs. The input SNR for this Figure is 0 dB, and it can be 
seen that the loop requires about two and one half minutes to lock to this low- 
levet signal. The loop requires longer times to lock to lower SNR signals and 
shorter times for high SNR inputs (e.g. about 15 minutes for -15 dB and 2 minutes 
for + dB SNR). 

Figure 13 shows the steady-state response of the SMAPLL2 to the desired 
range of SNR’s of step inputs. Figure 14 gives the same data for the 300 knot 
ramp inputs. The vertical axis is the standard deviation of the SMAPLL2 output 
estimates of phase in units of cec. Both of these Figures show that the loop s 
steady-state error will be less than plus or minus 3 cec 62.5% of the time whether 
airborne at 300 knots or on the ground motionless. 

A further demonstration of the improvement in signal variance provided by 
the SMAPLL2 is provided by plots of the histograms for signal plus noise input 
phase and the loop output estimates of phase. Figure 15 gives the input 
histogram for a -10 dB SNR case and Figure 16 gives the loop’s output histogram. 
The improvement made in signal variance by the 5MAPLL2 is obvious from this 
Figure, and is directly proportional to the SNR improvement through the loop. 

Ail SMAPLL2 parameters have been optimized for tracking in the noisy 
-20 dB SNR case and with the 300-knot velocity constraint imposed. Performance 
is considerably enhanced by higher SNR’s and slower velocities. The SMAPLL2 
routine requires about 150 bytes of microcomputer memory and takes about 200 
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Figure 12. SMAPLL2 Transient (Lock-Up Mode) Response to 180® Step 
(Top) and 300-Knot Ramp Inputs (Bottom) of 0 dB SNR, 
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Figure U Steady-State Response of SMAPLL2 to 300-Knot Ramp Inputs: Standard deviation 
of loop outputs in cec vs. range of input SNR's. 
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microseconds io eKecute. Appendix E gives the details of the simulation and 
additional plots describing the loop's performance. Appendix F gives the nticro- 
computer program listing for the SMAPLL2. 

In addition to the interrupt service routine, auto-sync subroutine, and 
SMAPLL2 subroutine already described (and the data output subroutine to be des- 
cribed in the next chapter), the OSP employs a “main", or idling, routine. At 
system turn-on time this routine initializes all necessary memory locations to their 
proper contents; it then simply idles (wasting time) until it is Interrupted by the 
microcomputer interface module with new Omega data. This main routine would 
take the form of an airborne navigation program in a full-blown ONS. In its 
current form it requires 100 bytes of memory to initialize OSP parameters and 
takes about 100 microseconds to execute. 

By accumulating the CPU time and program memory requirements for alt 
the OSP algorithms presented, it is possible to estimate the CPU time and program 
storage currently unused (that which is available for future navigation processing 
algorithms). After system initialization and automatic synchronization during the 
first 100 seconds after turn-on, the SMAPLL2 requires about 200 microseconds 
after each Omega interrupt every 10 milliseconds. This leaves about 98 /o of the 
available CPU time currently unused (i.e. available for navigation processing). 
AM OSP algorithms, including the Data Output Routine described in the next 
chapter, require o total of about 1000 bytes (1 K) of program memory. There- 
fore, about 3K of the available storage is currently unused and is available for 
the navigation algorithms. 
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III. OSP OUTPUTS AND USES 


A. Output Data, The Omega sensor processor receives the Omega signal 
off-the-air, performs automatic synchronization to the Omega transmission format, 
and tracks the phase of each receivable station. The data which the OSP makes 
available is the relative phase of each transmitting station and an indication of 
the relative signal strength of each station's signal with respect to the atmospheric 
noise level. This is the data needed by a navigation processor to derive the 
receiver's posi;ion« 

The microcomputer-based OSP makes maximum use of the inherent design 
flexibility afforded by the software by utilizing a data output subroutine. The 
outpu'- routine is flow-charted in Figure 17, and the actual program is listed in 
Appendix F. The routine is built around three user-inserted code words, as 
illustrated in Figure 18. (Appendix D explains the insertion of data into the 
microcomputer). One bit of the "option" code word selects a single output or 
two outputs to be time multiplexed to the output latch (this latch is described 
in the hardware section). Another bit selects output of the phase or signal 
cjuality of the channels chosen by the second and third code words. The last 
six bih of tne option code word dictate the multiplex time in tens of seconds to 
be spent on the output of each channel's data (if the multiplex option was selected 
by bit one). The second and third code words select the stations for which data 
is to be displayed. Each code word is treated as two four-bit station identifiers, 
and each station Identifier is a number between zero and eight (corresponding to 
zero and the Omega stations A through H). 
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MUX (Bit D^) is set = 1 for multiplexing tv*/o outputs. (D^ = 0 will result in 
LOP(l) or S/N(l) output only.) 


S/N or LOP (Bit D^) is set = 1 for LOP cutput (D^ = 0 will result in S/N output). 

MUX Time (Bits - Dq) is set to a binary number corresponding to the number of 
10 second intervals to be spent outputting each LOP or S/N before switching to other 
LOP or S/N. 


LOP and S/N Addresses 

Code Words (K)01 and 0002 


StG'ion I (Bits “■ D^) is o binary number between 0 and 8. The number 0 selects 
clock drift on!)/; while 1-8 select Omega stations A-H. 

tOP 13 genercied taking phase of station 1 - station 2. 

if S/N is selected by making bit of option code word » 0, then the S/N of 
station 2 is output. 

Figure 18, User- Inserted Code Words for Outfit Routine. 
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If fhe phase option has been selected, the numbers in the first channel 
code word represent the phase of station minus the phase of station to be 
output. If the number inserted for station ^'"‘l is zero, and the number of station 
^2 is 1-8, the phase of the stai-ion selected in station "2 with respect to the OSP 
TCXO will be output. This measurement of "single station phase" Is exactly that 
which has been tracked by the SMAPLL2 tracking filter described in the software 
section. An example of the single station phase output appears in Figures 19 
and 20. Figure 19 shows the phase of a strong SNR station (as received in Ohio) 
versus the OSP clock. The Omega signals are phase-locked to a stable atomic 
standard at the transmitter, so that the constant velocity offset in phase is caused 
by an offset in frequency of the OSP clock. Figure 20 shows the single station 
phase of a station whose received SNR is lower than in the previous graph. 

This signal also exhibits the ramping phase measurements, but the track is less 
constant or noiser than the strong station's plot. 

If the phase option has been selected and neither of the numbers in the 
first channel word are zero, the phase difference of station ^1 minus station 72 
will be output. The phase difference of two Omega signals gives a line-of-position 
(LOP) as explained In Appendix A. Figures 21 and 22 show the LOP's plotted for 
various station pairs as received in Ohio. The phase ramping noted in the single 
station phase plots is eliminated in the LOP's J>y virtue of the fact that subtracting 
two phases measured against the same clock reference effectively cancels the 
clock from the measurement. Figure 2la shows the strong station pair C-D, 
and also illustrates the predictable diurnal phase shifting effects caused by sunrise 
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Figure 19. Single Station Phase of North Dakota Minus Clock. 
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Figure 20. Single Station Phase of Trinidad Minus Clock 
and Argentina Minus Clock. 
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Figure 21a. 24 Hours Diurnal Data: Hawaii - North Dakota LOP. 
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Figure 21b. 24 Hours Diurnal Data: North Dakota - Trinidad LOP. 
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Figure 22. 24 Hours Diurnal Data; North Dakota - Argentina LOP. 
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and sunset. Figure21b shows the D-GLOP which is anoisier pair, but one which 
is still highly usable for navigation. Figure 22 shows the D-F LOP which is very 
noisy, and is mostuseful for navigation during the atmospherically "quieter" night- 
time hours (when the propagation path is in darkness). 

For means of comparing the software-based OSP with an all hardware 
OSP, Figure 23 is given which shows LOP output's for both designs. This chart 
illustrates that both OSP designs track the Om^a signal when in a ground-based 
environment. However, the second-order tracking loop employed in the software- 
based OSP enables the tracking of high speed aircraft with more SNR improvement 
than is possible with the first-order hardware version. This is verified by the computer 
simulation described in Appendix E. 

If signal quality is selected in the option word instead of phase, the output 
will be the signal quality of station ^2 of the first channel word. Figure 24 shows 
plots of the signal quality for several stations. Note that this gives only a relative 
SNR indication in its present form; although this measurement could be calibrated 
to give true SNR. 

The multiplex option represents a means for time-sharing a single chart 
recorder or other output device. If this option is selected in the option code word, 
the output routine will multiplex phase or signal quality (whichever is selected) 
of the first and second channel code words to the output latch. The multiplex 
time, or time spent outputting the data for each channel word, is determined 
from the last six bits in the option word. This six-bit binary number represents the 
number of ten second intervals to be spent on each channels data before switching 
to output the data indicated for the other channel word. Figure 25 gives an example 
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Figure 23. Omega LOP Output Traces -- Software Loop Output on Bottom. 
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Figure 25. Multiplexed LOP Outputs (90 Seconds MUX Time) 
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of a multiplexed LOP output. 

B. OSP Uses. The primary intended use for the OSP with the outputs 
as presented here is to provide Omega data to an airborne navigation processor. 

This navigation processor is then to derive pilot-usable information such as 
current position, course deviation from desired course, and miles-to-go to 
destination. For this use the OSP is thought to perform well. 

For the airborne use of OSP data a sophisticated navigation processor 
is needed; however, there are applications for the use of OSP data which requ -e 
little or no navigation processing. The boat navigator can observe the LOP out- 
puts from the OSP, apply a skywave correction available in published tables , 

[ 17 ] 

and find his current position on published Omega charts. In this slow-speed 

application, the boat navigator can perform his own navigation processing by 

calculating the course and distance-to-go from his charted. Omega position. . 

Many applications of Omega data utilize remote computer facilities to 

process the raw Omega phase data from an OSP. In these applications the Omega 

data is collected via an OSP on board an object to be tracked (e.g. weather 

balloon, radiosonde, water bouy, life raft, etc.). The Omega coordinates are 

then retransmitted via some other communications link to the remote tracking 

. . . r [18,19,20] 

station for processing and cojjversion to position information. 

The U.S. Coast Guard is the U.S. government agency responsible for 
maintaining the Omega network in navigational ly-precise, working order. 

This entails providing a network of Omega monitor stations to receive and record 
the Omega signals at several geographic locations around the world for transmitter 
fault detection purposes and signal propagation studies. The OSP presented herein 
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Is ideal for monitoring purposes such as thiso The data can be stored on a digital 
or analog output device for later collection, or it can be transmitted to the 
monitoring headquarters via phone lines in near real time. 

Finally, the sihgle station phase measurements such as Illustrated in Figure 
19 can be used to calibrate laboratory oscillators to near atomic clock precision. 
Since the Omega signals are locked to atomic standards at the transmitters, the 
received signal versus local clock measurement shows the offset of the local clock. 

By adjusting the frequency of the local clock until the phase ramping gives way 
to a straight line (constant phase measurement) the local oscillator can be calibrated. 

IV. SUMMARY AND CONCLUSIONS 

A microcomputer-based low-cost Omega sensor processor for general 
aviation use has been designed, implemented with minimum hardware and 
innovative software and tested. The sensor processor hardware configuration 
has been presented and the software has been described. 

A FORTRAN simulati analysis of the functions of automatic synchroni- 
zation and PLL phase tracking has been performed to optimize routines to be per- 
formed by the OSP microcomputeri A data base of real Omega data was collected 
using the system's receiver and interface modules and was used as the simulation 
input to obtain realistic results. 

The microcomputer routines developed have been shown to output data 
comparable to an all-hardware OSP and, in addition,the software— based OSP 
offers data options and capabilities not available with the all-hardware model. 

The OSP provides outputs suitable for use by a navigation processor, and the 
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system presented includes nearly all the hardware needed for an airborne Omega 
navigation system. 

The feasibility of Implementing a minimum hardware software-based OSP 
of higher performance and lower cost than is currently available has been demonstrated. 
Further research is encouraged to develop navigation processing microcomputer soft- 
ware to complete the Omega navigation system tor general aviation use. 

V. RECOMMENDATIONS FOR FUTURE WORK 

The OSP presented herein represents nearly all the hardware needed for a 
complete Omega navigation system. The software OSP routines provide Omega data 
usable by a navigation processor to derive pilot-usable information. Needed to 
complete the system for general aviation use is a microcomputerized area navigation 
routine. One possible navigation technique using Omega data has been functionally 

[213 

described by Zervos. In addition, a propagation correction routine or some 
other means Tor applying skywave corrections is needed. More work also needs 
to be done to develop a lower-cost, effective pilot display. 

A software-based navigation system is afforded inherent flexibility in that 
the computer memory can be reprogrammed to perform other functions. The OSP 
design presented can be used for some other navigation system by augmenting or 
replacing the Omega front-end nrfodule with a LORAN, GPS, or other front-end 
and modifying the software. The potential for a hybrid Omega-LORAN system 
seems great, considering the possible sharing of system hardware and software by 
the two systems. '| 
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Aside from the noted receiving system developmental work, additional research 
Is needed to determine means for improving the redundancy and reliability offered by 
the Omega transmitting network. Whether this takes the form of a supplemental 
navigating system in the receiver-computer or an actual transmitter station aug- 
mentation, pilot confidence in on area navigation system must be justified and 
earned. Thereafter, a means for efficiently integrating Omega into the National 
Airspace System needs to be pursued. 

VI. ACKNOWLEDGEMENTS 

The author wishes to thank Mr. Jack Reid of the NASA/Langley Research 
Center for his support of this project as the Project Monitor of the NASA Tri- 
University Program. Gratitude is expressed to Dr. Richard McFarland, Director 
of Avionics at Ohio University for his support and encouragement. Dr. J. E. 

Essman of the Electrical Engineering faculty is also acknowledged for his con- 
structive advice. 

Special appreciation is reserved for Mr. Ralph Burhans, Avionics Research 
Engineer, and Dr. Robert Li I ley. Assistant Director of Avionics, without whose 
technical expertise and enthusiastic support this thesis would not have been possible, 
Mr. Kent Chamberlin and Mr. Paul Blasche, Graduate Research Associates, have 
the author's sincere thanks for their technical consultation. 

The author's Tri- University Program student colleagues at Ohio University, 
Princeton University, and MIT are also acknowledged as being instrumental in 
this project's completion. Thanks also go to Mrs. Roma Beverage and Mrs. Hope 
Mills for their work on the manuscript. 


- 49 - 



VII. 


REFERENCES 


Cl] Department of Transportai-lon, U.S. Coast Guards Omega Navigation 
System User Manual^, CG-ONSOD^ 2100 2nd St. SW, Washington, 

D. C., preliminary tlraft in January, 1976, 

C23 Pierce, J.A. , et.al, . Omega; A World-Wide Navigation System; 

System Specification and Implementation , AD 630 900, National 
Technical Information Service, Springfield, Virginia, issued May, 

1966, 

C33 Morris, Peter, and Milton Cha, Omega Propagation Corrections; Back- 
ground and Computational Algorithm , Report No, ONSOD-01-74, 
Department of Transportation, U,S, Coast Guard, Washington, D. C. , 
December, 1974. 

[43 Radio Navigation Group, Northrop Corporation, Omega Navigation 
Set AN/ARN-99 (XN-2), Phase II, Final Engineering Report , 2301 
W. 120th St., Hawthorne, California, December, 1973. 

C53 Ohio University Avionics Engineering Center, Prototype Omega Receiver 

for Use in Data Collection and Evaluation , Athens, Ohio, September; 1976, 

[63 Burhans, R.W., Low-Cost, High-Performance, VLF Receiver Front-End , 
Technical Memorandum (NASA) 18, Avionics Engineering Center, ~ 

Ohio University, Athens, Ohio, September, 1975. 

[73 Wright, Lee, Improvements for Omega RF Preamplifiers , Technical 

Memorand'.'m (NASA) 24, Avionics Engineering C^teT, Ohio Univer- 
sity, Athens, Ohio, April, 1976. 

[83 Lilley, R. W., and R. J. Salter, Simulation Analysis of a Microcomputer- 
Based, Low-Cost Omega Navigation System, Proceedings of the Bi- 
centennial National Aerospace Symposium of the Institute of Navigation, 
Warminster, Pennsylvania, April, 1976. 

[93 Lilley, R. W. , and R. J. Salter, A Microcomputer- Based Low-Cost Omega 
Navigation System, Prpceedings of the First Annual Meeting of the Inter- 
national Omega Association, Arlington, Virginia, July, 1976. 

[103 Salter, R. J. , “Navigation with Mini-O, Part 3, Software”, Byte Magazine , 
Peterborough, New Hampshire, \|lol. 2, No, 4, April 1977. 

[113 MOS Technology, Inc., MCS65Q0 Microcomputer Family Hardware Manual , 
950 Rittenhouse Road, Norristown, Pennsylvania, August, 1975^ 


- 50 - 



1 


C12] MOS Technology, InCo, MCS6500 Microcomput-er Family Programming 

Manual, 950 Riti-enhouse Road, Norristown, Pennsylvania, August, iy75» 

[13] Chamberlin, K. A., Digital Correlation Detector for Low-Cost Omega 
Navigation, Technical Memorandum (NASA) 19, Avionics Engineering 
Center, Ohio University, Athens, Ohio, February, 1976. 

[14] Raab, Frederick H,, and Jerome R. Waechter, The Use of a Counting 
Phas e Detector to Preprocess VLF Radionavigation Signals, to be 
published in IEEE Transactions on Aerospace and Electronics Systems. 

[15] Chamberlin, K. A., and R. W. Li I ley. The Memory- Aided Phase- 
Locked Loop: A Model and Proposed LSI Chip Implementation , 

Final Report prepared for Naval Avionics Facility, Indianapolis, 

Indiana, No. EER26-1, Avionics Engineering Center, Ohio 
University, Athens, Ohio, March 1976. 

[16] H.O. Pub. No. 224 (111 -C) Series: Omega Propagation Correction 
Tables, U.S. Naval Oceanographic Office, Washington, D.C., 1971. 

[173 Omega Navigational Charts, 7500 Series , U.S. Naval Oceanographic 
Office, Washington, D.C. 

[18] Beukers, John M., Windfinding Using Navigational Aids , presented at 
The Third Symposium on Meteorological Observations and Instrumentation 
of The American Meteorological Society, February, 1975. 

[19] Baker, Donald L., and Chfistopher S. Welch, A Prototype Legrangian 
Current Buoy Using the Carrier Plus Sideband (CSB) Retransmission of 
Omega Navig^ation Signals , Proceedings of the First Annual Meeting 
of the International Omega Association, Arlington, Virginia, July, 

1976. 

[203 Rupp, Walter E. , GRAN: Experiences and Activities , Proceedings of 
the First Annual Meeting of The International Omega Association, 
Arlington, Virginia, July, 1976. 

[21] Zervos. Arthouros, K., A Performance Simulation of a Potential Class 
of Low-Cost Omegg/Air Data Navigation Systems, Masters Thesis, 
Princeton University, Department of Aerospace and Mechanical 
Sciences, September, 1975. 

[22] JOLT CPU Hardware Manual a nd JOLT DEMON Software Manual, 
Microcomputer Associates, Inc., Santa Clara, California, First 
Edition, 1975. 


- 51 - 


123] Lilley, Robert W., An Assembler for the MOS Technology 6502 Micro- 
processor QS Implemented in the JOLT (TM) and KIM-1 (TM)/ Technical 
Memorandum (NASA) 44, Avionics Engineering Center, Ohio Univer- 
sity^Athens, Ohio, November, 1976. 

1243 Schwartz, Mischa, information Transmission, Modulation and Noise, 
Second Edition, McGrow Hill Book Company, New York, 1970, 



VIII. APPENDICES 


A, The Worldwide Omega Navigoi'ion System . The worldwide Omega 
System is a network of very low-frequency (VLF) transmitters located at eight 
locations around the world/ as shown in Figure A-1 . Each station transmits for 
approximately one second out of each ten seconds on the frequencies of 10.2/ 

13,6, and 11.333 KHz. As shown in Figure A-2, a unique transmission pattern 
is formed hy the scheduled length of each station s transmission. 

Omega was originally conceived by Professor J. Pierce of Harvard Univer- 
sity following his successful demonstration of the long-range phase stability of 
VLF signals in the 1940's. Originally developed by the U.S. Navy for use on 
seagoing ships. Omega has only recently been used in airborne applications. 
Operational responsibility for the system was transferred to the U.S, Coast Guard 
in 1971, a move indicative of increasing civilian use of Omega, As of this writing, 
seven of the eight permanent transmitters are in operation, and the Coast Guard 

I 

is beginning implementation of approximately forty monitor receiving stations 
around the world. 

Since the wavelength of a 10.2 KHz signal is about 16 nautical miles, 

.an effective antenna for radiating this frequency must be on the order of several 
miles long. Two types of antenna structures are utilized for Omega transmission. 

As diagrammed in Figure A-3, these are the tower and valley span antennas. 

Signals in the VLF region propagate in the waveguide mode using the conducting 
earth and ionosphere as waveguide bounds. As shown in Figure A-4, a single 
waveguide mode dominates when several hundred miles from a transmitter, and 
due to this single, stable mode of propagation the VLF signals can be received in 
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Figure A-2, Omega Signal Transmission Format 


Figure A-1 . Omega Station Locations 
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Figure A-3. Omega Transmitting Antenna Structures 
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excess of 6,000 miles and at all altitudes (even beSween mountains and beneath 
the surface of water). 

The height and smoothness of the ionosphere boundary is changed by vary- 
ing amounts of radiation received from the sun, and this effective change in wave- 
guide height results in a signal phase velocity change. The phase velocity change 
occurs when the signal propagation path to the receiver changes from light to 
darkness and vice-versa, and is, therefore, referred to as "diurnal" shift. This 
diurnal shift is highly predictable, and the received phase can be adjusted by 
applying a skywave correction (SWC) factor which is unique for month, day, 
time-of-day, and receiver location. Another highly predictable factor is the 
varying conductivity of the earth waveguide boundary. Ground conductivity 
maps have been prepared and a propagation correction can be applied for particular 
receiver location. Ionospheric disturbances caused by sunspot activity and solar 
flares are not as highly predictable and limit the position-fix accuracy attainable 
with Omega. 

Techniques for the reception of the Omega carrier signal differ; however, 

detection of the signal's phase Involves comparing the received signal to a locally 

generated signal of the same frequency. Differing degrees of sophistication are 

possible in the phase detection process; more sophisticated meaning more accurate 

12 , 

local clock signals. If a clock frequency on the order of a part in 10 (atomic 
standard accuracy) is available, and it is possible to synchronize the phase of 
the local clock to the time standards employed at the Omega transmitters, the 
Omega system can be used in the rho-rho positioning mode by simply measuring 
the tlme-of-arrivai of two or more Omega signals (propagation time being 
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proport'iona i to distance traveled)® If a local clocU of this high accuracy is avail- 
able# but it is not possible to synchronize iis phase with the Omega time standards, 
it is still possible to navigate in the rho-rho mode® However, it is first necessary 
to determine the receiver's position within height nautical miles (since a cycle at 

10.2 KHz is 16 nautical miles long and each phase repeats every 16 miles). 

This is known as lane resolution, and several techniques exist for resolving lanes 
to within needed limita. 

If a local clock of atomic standard accuracy is not avi^Mable, a less stable 
(and less expensive) clock may be used provided short-term stability constraints 
are obeyed. For example, a clock stabi I iSy of one part in 10 per day will drift 
no more than one centicycle (1/100 of a 360 degree cycle) in ten seconds. There- 
fore, if this local reference clock is, in effect, used within each ten second Omega 
frame separately, very little accuracy is sacrificed. Specifically, phase differenc- 
ing can be performed by first detecting each stat!on°s signal with respect to the local 
clock and then subtracting two phase measurements. This effectively subtracts 
the clock from the measurement. Each phase difference is a point on a hyperbolic 
contour. Geometrically, lines of constant phase difference form hyperbolic contours 
referred to as hyperbolic lines-of-position (LOP“s), as shown in Figure A-5. 

Navigation by reference to these LOPE'S is known as the hyperbolic mode, end it 
is again necessary to determine Initial position to within one-half lane (lane resolution). 
It is necessary to use iwo LOP°s (at least three stations) to determine position, and 
greatest accuracy is achievable when using LOP“s which Intersect at close to a 90 
degree angle. 
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Paramounf In the detection and estimation of Omega phase is the problem 
of noise. Amplitude noise in the VLF region is caused primarily by lightning and 
is made up of both gaussian distributed noise (a conglomeration of many distant 
lightning sources of various strengths and disiances) and impulsive noise (closer 
lightning sources). Transmitted power of all Omega stations Is ten kilowatts^, 
and a representative sample of received signal and noise amplitudes appears in 
Figure A-6. Also of importance is the noise caused by precipitation upon the air- 
craft and antenna. These additive amplitude noises cause the detected phase 
to be contaminated with noise (i.e. phase jitter). Due to the several contributing 
noise sources and allowing for less than perfect propagation corrections, the stated 
Omega system position fix accuracy is + - one nautical mile (daytime) to + - two 
nautical miles (nighttime). This basic accuracy limitation constrains Omega use 

to an area (non-precision) navigation aid. 

B, Omega Data Base. For purposes of obtaining realistic data for a 

simulation analysis of Omega sensor processor (OSP) functions, a data base of real 
Omega ground and flight data was collected. So that the data used as input to 
the simulation would be the same as that seen by the OSP microcomputer algorithms, 
the data was collected using the OSP hardware front-end and interface modules 
which had been fabricated and perfected before the OSP software was developed. 

For purposes of evaluating the performance pf the simulated automatic 
synchronization and PLL tracking processes it was necessary to collect Omega 
data under the following conditions: (1) ground-based receiver with noise conditions 
found throughout a 24-hour day; (2) ground-based receiver with and without the 
overwhelming strong station (North Dakota) being received; and (3) airborne 
receiver in a general aviation aircraft. 
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The Omega phase data was collected via the equipment set-up as shown 
in Figure B-1. The microcomputer provides an interrupt signal when new data 
is latched and ready to be read* For data collection the 100 Hz interrupt signal 
was used as a magnetic tape recorder "write" signal (one phase sample written 
every ten milliseconds). In addition, it was possible to collect data only during 
every other ten second time frame. This was to allow time to writ3 an end-of- 
record on the tape to separate the ten second data frames. Note that this is 
wideband 10.2 KHz Omega phase data taken directly from the 30 Hz bandwidth 
receiver front-end and sampled by the interface module every ten milliseconds. 

The ground-based data was collected at Ohio University°s Clippinger 
Laboratories, A 6,1 meters long wire^ mounted on the roof of Clippinger Labora- 
tories was used a", the receiving antenna. An active preamplifier was colocated 
at the antenna and feeds the signal through approxim iteiy 30,5 meters of coaxial 
lead-in cable to the receiver in the Omega Laboratory, The electrical environ- 
ment of this receiving antenna is one of considerable interfering noise. The noise 
is primarily from two sources; (1) 60 Hz harmonics radiated by high-powered 
electrical equipment in the building; and (2) audio modulation of a nearby AM 
radio transmitter. Therefore, the ground-based data is of worst-case variety. 

The data without strong station North Dakota was collected as that station was 
going off-the-air for scheduled maintenance. 

The airborne data was collected with the equipment set-up of Figure B-1 
installed onboard the Avionics Engineering Center DC-3 flying laboratory. The 
data was collected outbound to Chicago O' Hare International Airport from Ohio 
University Airport at Albany, Ohio (and also inbound on the return flight). The 
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Figure B-1. Omega Data Collection Setup, 





outbound flight was conducted between the hours of 7:00 and 8:30 AM on October 20, 
1976, The cruise altitude was 2438 meters and VFR conditions persisted through- 
out the flight. The return flight was between the hours of 6:00 and 7:00 PM on 
the same day and the cruise altitude was 2743 meters in VFR conditions. The 
aircraft's ADF sense antenna mounted on top of the fuselage was used with the 
preamplifier mounted on the ceiling inside the aircraft. The preamplifier fed 
approximately 6.1 meters of coaxial lead-in cable to the receiver. 

All data was recorded on magnetic tape in the binary format presented 
by the microcomputer interface module. The phase of the Omega signal is 
recorded every ten milliseconds, each data point being quantized into the 
interval from zero to 255. The OSP“s temperature compensated crystal oscillator 
(TCXO) was used as the local reference clock; therefore, the slight frequency 
offset of the TCXO appears as a constant-velocity phase ramp in the data. 

Figure B-2 summarizes the contents of the data base. 


File ^ 

Type 

Data 

Dote 

Collected 

Location 

Stations 

Transmitting 

1 

Ground (24 hours) 

January, 1976 

Athens, Ohio 

B,C, D, G 

2 

Ground (D off air) 

June, 1976 

Athens, Ohio 

B,C,G 

3 

Ground (D off air) 

June, 1976 

Athens, Ohio 

B,C,G 

4 

Flight (VFR) 

October, 1976 

Athens-Chicago 

B,C,D,G 

5 

Flight (VFR) 

October, 1 976 

Chicago- Athens 

B,C,D,G 


Figure B-2. Omega Data Bose Summary. 


c. Omega Sensor Processor Circuitry. The Omega sensor processc, is 


composed of three distinct hardware modules and a general-purpose microcomputer. 
The circuit details of the Omega preamplifier, receiver front-end, and micro- 
computer interface modules are given in Figures C-1 , C-2 and C-3, respectively. 
The microcomputer system is described in Appendix D. 

D. JOLT Microcomputer System . 

1. Hardware. The JOLT Microcomputer System is available from 
Microcomputer Associates, Inc., Santa Clara, California, As implemented in 
the Omega Sensor Processor, the JOLT system consists of a CPU board and a memory 
board. The CPU board containing the 6502 CPU chip, 6530 chip (PIA, TTY 
Monitor ROM, and Interval Timer), 6520 PIA chip, TTY interface circuitry, 
and 512 bytes RAM is available in kit form for $150. The 4K RAM memory board 

is available in kit form for $200, Circuit diagrams for the CPU and memory boards 

C223 

appear in Figures D-1 and D-2. 

As can be seen from these diagrams, the JOLT incorporates a 16-line 
address bus and an 8-iine data bus, along with an adjustable 700 KHz - 1 MHz 
oscillator. System architecture includes one accumulator, two index registers, 
one processor status register, one stack pointer register, and one program counter 
register. All are 8-bit registers except for the program counter which is 16-bits 
long. 

The system accommodates two fully vectored interrupts. One is a non- 
maskable interrupt (NMI), and the other Is an interrupt request (IRQ) which can 
be disabled and thereafter ignored by the software. 
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Figure C-1. Dual Purpose Preamplifier Module. 
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Figure C-2a. Receiver Front-End Module, Two-Stage 
VLF Bandpass Amplifier. 
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Figure C-3. Microcomputer Interface Module Circuit Diagram 













Figure D-1. JOLT CPU Board Schematic Diagram, 









Figure D-2. JOLT 4K RAM Memory Board Schematic Diagram. 
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The JOLT system requires approximately 2 amps at + 5 volts and 50 milli- 
amps at -10 volts DC. After these supplies have been turned on, the user need 
only push the system reset button and the system is ready to accept user-commands 
from a TTY keyboard. 

2. Software. Figure D-3 gives the instruction set for the JOLT 
system and details about each instruction. It can be seen that the average in- 
struction fetch and execute time is about 3 microseconds. The 52 instructions 
combined with their various addressing modes provide one hundred and fifty unique 
operations. Figure D-4 gives more software details and Figure D-5 gives the 

*^223 

JOLT system memory map. 

The TTY monitor routine which is resident on the 6530 IK ROM chip 
allows the user to operate the JOLT from a TTY keyboard. Commands include! 

(1) display contents of memory^ (2) display contents of Cl* registers; (3) alter 
contents of memory or CPU registers; (4) begin program execution; and (5) punch 
and read paper tape. In addition, several monitor subroutines can be called to 
type data on the TTY printer. 

In addition, a 6502 Assembler vyhich can be executed via a larger computer 

[ 23 ] 

system employing FORTRAN IV, has been obtained. This Assembler capability 
allows the programmer to write microcomputer programs in mnemonics rather 

than hexadecimal number codes, and allows for use of variable names and labels 

rather than numerical codes and locations. 

E. FORTR.^N Simulation Analysis . In order to test the feasibility of per- 
forming OSP functions of gutornatic synchronization and PLL phase fracking in 


-71- 







* (■ClUVwlO* 

/ yco>r>tO 

«tOt ««OD«*‘tO 
•t, MMO«*t)r} 


A ACCUMUVAtO* 

M «M*( 04 y nn AOOMSS A 

Mk MtMOOV Af« <t*c« »»NTt« v_ 


MtTRACt 

AnO 


•K rAA6t i 


[ieAA»i.V] 


r AWot uTi wcnKfP jHti'»rctrvtAOO»i:« » 5 »oi«no*TAOOH.GT»»ii*«f« 
VjfTHi’ilCONh iwo th:» 0 •*»! 0» r**t 

HflfKlQiMOtAeC? TH< ^COHC iTTl 0» T*** rtrttRliCTKJN.S AOOIO TO 
(MtOCK D»SC*l»I ^*»0 TMt CARRY Tm| «fV--lT% AOtHTl TO A lOCATlON CN 
►ACtZlIWWHiCNCOATAiMST*^ •lOWO«'J*«rHO* *»l lA TMfN|RT«VT| 


THC Off RAMO IS CONt R>*HO M TMl UCOMO RYTt Of 


ADO»W> «T»Q - T»*C UCOMO pyTt Cf T»*t i*i|TRuCTiO»ieOHTA««S tHt 
RlORlOROfRima* THI IfMCrrvt AOORIl' TKCTMiROBrTtCO«(TA.»rtTM* 

• HIGH ORMR BITS Of ThI EffCCTiVt AOORiSC 


I ffRQ * 00 «HO» *C - MCONOtTtl tO»iTAtfiiTHtBU)*ORO€Bi*TSOf t* 4 f 
CffCCrry^ AOORCSS Th( BhiChOROCII nn ARE ZERO 


ORE BYTE IRSTRuCTKMOfERATIRGOft the ACCiAIuLATOR 


AOOrOTOtM(ii«OC^I XARRYiSOROftSOi TO rORW THE ION OROfR BYTt Of 
TMCIA THSMICMOROfRBYTCOf TMCEAttZfROf 


Setter MdS Te^nology 6502 CPU 


Instruction 


r 

1 * 0 ^ ^ 
Jaa 


c 

; *, 

1 " 


> 





-. 1 ww 
^yjSa: 



• 

ismi 

R>l 






• 


L*i‘ 

1 . 

1 c 

■ 




COMMAND 

R 

M 

a 

• 

G 

WH 

LH 


OPERATION 

Print out the CPU registers. 

Print out contents of memory. 

Alter contents of CPU registers or memory. 
Begin program execution. 

Write contents of memory in hexodecimctl. 
Read punched paper tape. 


Figure D-4, JOLT System TTY Monitor Commands Summary 
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software, a FORTRAN stmulaHon analysis was performed. This simu lotion also 
served to optimize parameters in both algorithms before they were assembled in 
microcomputer code. The input data for much of the simulation analysis was 
the real Omega data base described in Appendix B; however, since the specifi- 
cation of system performance requires complete and absolute knowledge of the 
input data, computer generated data was also used. 

Each simulation is summarily described, and the actual listing for each 
FORTRAN simulation appears at the end of this Appendix. 

1 . Noise Modeling. Since the real Omega data base contained 
uncertainties such us signal-to-noise ratios (SNR), phase measurements designed 
to approximate the data at the output of the OSP front-end and microcomputer 
interface modules were generated with the computer. The Omega noise is both 
gaussian and impulsive. Gaussian noise is easily generated by the computer and 
was judged to sufficiently approximate real Omega noise following comparisons 
of statistics of the real and modeled noise. 

The rtoise model must simulate the gaussian distributed input noise being 
bandpassed through the 30 Hz narrowband Omega receiver front-end module. 

The basis for such a model of a sinusoidal signal embedded in gaussian noise and 
received through a narrowband receiver is developed by Schwartz. 

The signaT developed across the antenna can be expressed as: s(t) = 
(x+A)cos(Art'l'o)+(y)sin(wt), This expression neglects possible interferring carrier 
terms from other Omega frequencies and other strong signals near the desired 
frequency. Acos(wt+o) is the desired signal term, and x and y are the white, 

2 

gaussian noise terms (with x and y independent, of zero mean, and variance o ). 
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These noise terms can be generated by calls to system library random number 
generator routines, RANDNR and RANDNRB. The narrowband receiver front- 
end module is simulated by a low-pass filter (LPF) subroutine, GRANGE. Pass- 
ing the random noise terms through the LPF results in correlated noise samples. 
The phase angle detected at the output of the LPF front-end can be simulated 
by taking 0 = arctan ( ■ The density function of phase obtained in this 


manner is; 


-s2 1 , 2p. 

f(9) = + iv — cos 0 e 1 1 + erf(s cos 0) ] 

' ' 2tt 2 ir 


2 . • a2 

where s is the power signal-to-noise ratio equal to — 7 ^ . 

2 °^ 

Comparison of the histograms and autocorrelation plots in Figures E-1, 

E-2 versus E-3, E-4 gives an example of how well the modeled noise approximates 


the real Omega noise. Lack of autocorrelation of the simulated narrowband noise 
is due to the fact that the LPF subroutine GRANGE simulates a single-pole RC 
low-pass rather than the multiple-pole low-pass filter provided by the real Omega 
front-end. However, the amplitude of the R(0) term (noise variance) closely 
approximates that of the real data, and the histograms are fairly well matched 
(remember that the signal-to-noise ratios for the real data are only approximations). 

All plots were generated via calls to the system library routine MMPLOTT. 

2. Automatic Synchronization Simulation. The automatic synchronization 

simulation is flow charted in Figure E-5. 

The real Omega data base was used as the input data for this simulation, 
and the subroutine DATA is first cal led to Insure a random starting point in the 
data base (simulating a random OSP turn-on time). After a completed 10-second 
frame of data is read, the data is ‘'mirrored". Figure E-6 shows the 10-second 
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Figure E-4» Aofocorrelotion of Real Omega Noise Phase for +20 dB SNRp 






























frame of daia and lis corresponding mirrored values plowed 10 points per line. 

Next, sample-to~sample differences are computed (see Figure E-7)« Each 
difference value is compared to the predetermined threshold constant, and a 
vector of "1” and "0" is generated which represents the correlated and uncorrelated 
samples. 

Before generating the vector of "1“ and ”0" by comparing phase dlffisrences 
to the threshold constant, it is helpful to smooth or average these differences. 
Several types of smoothing fi lter were tried, including straight averager, one 
and two-sided exponential, and median three, with various amounts of points 
tried in each filter. The most successful was found to be the 9-point straight 
averager. Using this filter, the accumulated sum of nine differences is compared 
to nine times the threshold constant. 

The threshold constant is most critical and has been determined (approxi- 
mately) using the decision theory approach as well as by simulation. The constants 

% 

determined by both techniques are in close agreement, thereby lending credence 
to the value chosen. 

Having determined this constant, the vector of "1" and “0“ is created by 
comparing each phase difference to it. Then the vector represents the received 
correlation pattern. This received correlation pattern can be compared (or cross- 
correlated) ”bit-by-bit‘' in a bilevel correlation process with the stored unique 

j ' ' 

Omega transmission pattern, A correlation counter is incremented for each 
coincidence of the stored and received patterns. Then the received pattern is 
shifted one bit and the entire cross-correlation process is repeated until all possible 
starting points have been tested. The number Of shifts necessary to yield the highest 



TOO cee 



Phase 

Differences 


0 cec 



Time « 10 sec/line 


Figure E-7, Sample-to-Sample Phase Differences Plotted Ten Points 
Per Line. 
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correlaHon counter value is assumed to be the desired starting point for the "A” 
time-slot of the Omega data frame. After shifting the original data the determined 
number shifts^, Figure 9 (in tejct) shovws the Omega phase samples plotted ten points 
per line. 

The routine's performance is improved if the operator inserts the stations 
which are expected to give coherent phase data samples (i.e. those stations 
received at or above 0 dB SNR). There was no noticeable difference in the routines 
ability to sync on ground or flight data, nor on data with the strong station North 
Dakota off-the-oif^ However, it is obvious that the technique depends on at least 
two stations being received above 0 dB SNR. Therefore, if only one station was 
sufficiently received the operator would have to identify which station it was. 

Note that if a two or three frequency system was aval table, only one station 
would be needed for sync (due to the unique lengths of adjacent time-slots). 

3. Phase Locked loop Simulation. Since the phase-tracking structure 
presented in the text (SMAPLL2) does not readily fit conventional models analyzed 
in the literature, computer simulation has been chosen rather than theoretical 
analysis to optimize loop parameters for the airborne Omega application. Computer- 
generated data has been used to optimize the parameters ond specify SMAPLL2 
performance. Finally, the real Omega data base was used to confirm expected 
tracking performance. Worst-case conditions of -20 dB SNR in the 30 Hz pre- 
detection bandwidth and 300 knots velocity of input signal phase in a direction 
radially to/from Omega stations have been simulated. 

There are two SMAPLL2 response modes of interest: transient (including 
lock-up mode) and steady-state response. The transient response is evaluated by 
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beginning with the loop in lock with the input signal; then step and ramp input 
signals are applied. The response to these inputs is observed for the entire range 
of -20 dB to +20 dB input SNR. Steady-stota ro’^ponse is evaluated by initializing 
the loop to a locked condition for various velocities of ramping input signals and 
calculating the loop's output error for each input signal after several hundred 
time-slots have been simulated. This response is also observed for the entire 
range of input SNR's, 

A three-step approach has been utilized to "home-in" on the best choices 
for the parameters of first-order bandwidth and second-order gain via this simulation 

(1) With the loop in the first-order configuration (SMAPLLl), 
the steady-state response for a range of first-order filter constants 
is examined. The value of loop filtering which results in an output 
standard deviation of or near zero over the desired operating range 
of input SNR is chosen. 

(2) With the loop in the second-order configuration {SMAPLL2), 
and using the chosen value of first-order filtering, the transient 
response is examined for a wide range of second-order gains. It 
Is desirable to choose the value of second-order gain which results 
in best transient response (fastest lock-up time). If this minimum 
lock-up time is excessive, the data rate from the first-order filter 
is too slow, and step must be repeated to choose a lower value 

of first-order filtering (i.e. to increase the loop's data rate). 

(3) With the chosen values of first-order filtering and second-order 
gain, the steady-state response of the SMAPLL2 to various velocity 
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rarnp inputs is examined, if the output error standard deviation is 
not within acceptable limits over the SNR operating range, steps 
’^1 and are repeated. A greater amount of first-order filtering 
is chosen while sacrificing some lock-up time. This trade-off is 
the design optimization that must be made to achieve the best 
combination of transient and steady-state performance for each 
particular application. 

Executing step of this approiich, the SMAPLLl yields the graph in 
Figure E-8. This graph gives the SMAPLLl steady-state response to zero-=velociSy 
input signals over the desired tracking range of SNR°s. Since a first-order filter 
of 5 bits (32 counts) (or .02 Hz bandwidth) is seen to give an output standard 
deviation of near zero cec over the operating range of SNR, this value was 
initially chosen. However, it was found that this amount of first-order filtering 
reduced the date rate to below acceptable limits (i.e. in step ^2 the lock-up time 
was found to be prohibitive). Therefore, the first-order filter constant of 4 bits 
(16 counts or .05 Hz bandwidth) was chosen. 

The .05 Hz first-order loop bandwidth yields acceptable Icck-up times for 
several values of second-order gain (again, determined by executing step '^2). 
Therefore, the best value of second-order gain is chosen as that which provides 
best steady-state performance for ramp inputs over the desired range of input SNR's 
(step ^'3)o Figure E-9 shows the steady-state response of SMAPLL2 to a zero-velocity 
input for several values of second-order gain. Figure E-10 gives the same data for 
a 300-knot input rdmp. It can be seen that a value of .25 for the second-order 
gain gives the best response for both input velocities. 
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Figure E-8. SMAPLL1 Steady-State Response to O-Velocity Input for First-Order Filter Equal to 
2 Counts (1), 4 counts (2)^ 8 counts (3), 16 counts (4) and 32 counts (5), 
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Figure E-9. SMAPLL2 Steady-St-ate Response to Zero-Velocity Input for Second-Order Gain Equal to 
Zero (1), .125 (2), .250 (3), .375 (4), .500 (5). 
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Figure E-9. SMAPLL2 Steady-StaSe Response to Zero-Velocity Input for Second-Order Gain Equal to 
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Figure E-IO, SAAAPL12 Steady-Stare Response to 30O^Knot Romp Input for Second-Order Gain 
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The transient response of the optimized SMAPLL2 to 180 degree step and 
300-knot ramp inputs is plotted in Figure E-M and Since a general aviation 


aircraft is usually on the ground for at least ten minutes after aircraft power is 
turned on, a lock-up time of this length is considered acceptable. 

The FORTRAN language listings for computer simu lations of Omego noise, 
Hie outomatlc synchronization algorithm, and the performance simulation of the 
SMAPLL2 are given as Figures E-13, E-14 and E-^15, respectively. 





l-]2. SMAPLL2 Trdhsienf (Lock*^Up Mode) Response to 180* Step(TGp) piid 
300-Knot Romp (Bottom) Inputs at +10 dB SNR. 
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C 

c 

e 

c 

c 

c 

e 
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c 

c 


G 

G 


ROUTINE COMPUTES THE r.UTOCORHEL&TION E(T) OF 30TH REAL OMEGA 
NOISE AND SIMULATED NOISE, AND PLOTS P.(T) ON SEPARATE GRAPHS. 

ROUTINE COMPUTES AND PLOTS THE HISTOGRAMS OF BOTH REAL OMEGA 
NOISE AND SIMULATED NOISE PHASE ON SEPAPATE GRAPHS* 

OPERATOR MUST SELECT SNR FOE SIMULATED NOISE 3Sf SETTING A. 
OPERATOR MUST SELECT SNR FOR REAL OMEGA NOISE BY READING OVEP 
DATA ON DATA BASE WHICH IS NOT WANTED. 

INTEGER TME 

REAL X{200) ,XX (10G) ,R (10G) ,XXA (100) 

Y(200) ,YY(20G) ,XXL(12) ,YYL(12) ,INTVl. . , 

XXL/'T’ , *1* ,* M' » ’E* ,* *,’M*, *S*, *E*,'C*/ 

YYL/'A», *a*,»T’,*0*»*C»,*0*,'R*,'3’,*U*,’L*,*^i*»*TV 

RANDOM NUMBER GENERATOR AND LPF TO GET SIMULATED NOISE, 
PARAMETERS ON LPF TO SIMULATE lOHZ BW AND IGOHZ SAMPLES, 


REAL 

DATA 

DATA 

GALL 

SET 


CALL 
CALL 
TAKE 
A: IS 
DO 5 


G 

C 


CALL INTGEN(5**11) 

DW-= 15. 

SMALLA=BW*12,57 
CAPA=2. «‘SQRT (4,/SM ALLA) 

INTVL=.01 
TME-2Q0 

CALL GRANGE (YY,INTVL,CAPA,S HALLA,TME) 

INTGEN /5**17) 

GRANGFtY,niTVL,CAPA,SMALLA,TME) 

ARGfAN OF RANDOM NUMBER RATIOS TO GET NOISE PHASE. 

THE SIGNAL VOLTAGE AMPLITUDE; POWER SNR IS (A*A) / (2*N*N) 
555 DO 5 1=1,200 
A=0.G01 

X (I) = ATAN2(YY (I) , {Y {i)*A)) 

CHANGE UNITS FROM RADIANS TO GEC, 

X (I) =X (I) /0.G628313S3 
5 CONTINUE 

GALL AUTGOR SUBROUTINE TO eOMPUTE AND PLOT R(T), 

GALL AUTGOR {X| 

PRINT 102 

FORMAT (IX, 'THIS IS THE 
SNR^1O.* ALOG10 (A*A*.5) 

BW2=2.*BW 

PRINT 101,SNR,3M2,INTVL 

FORMAT {IX, 10X,*SNR = ;,F5.2,* DB’,10X»’B« = *,?3.0,* HZ ,10X, 

C'SAMELiNG TIME = * ,F5.3,* SEC*) 

GALL TO HISTO SUBROUTINE PLOTS HISTOGRAM OF SIMULATED PHASF 

noise vs, phase IK GEC. 

CALL HISTO (X) , 

PRINT 10 5 

105 FORMAT (1X,» THIS IS THE HISTOGRAM OF SIMULATED NOISE PHASE VS. 

CPHASB IN GEC, •) 

PRINT 101 ,SNR,B»2,INTVL 




© ►d 
C 


102 


101 


AUTOCORRELATION FOR SIMULATED NOISE PHASE •) 


C 

C: 

c 


gRLI. MIA SUBRW PETORN REAL OMSGA PHASE PLUS NOISE* 

FIRST SiLFCT SN^^^^^ BY READING OVER UNWANTED DATA, 

500 COUTIMUEv,i;;;.v\-^-- ^ .-A 

DO 1 1=1,130 


Figure E-13. FORTRAN Program RSRAND Used to Analyze Real and Simulated 
Omega Noise. 
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CALL DATA (IOUT,S100) 

ISK.IP=I 

1 CONTI MUE 
DO 2 1=1,100 
CALL DATA(IOUT,6100) 

X(I) =FLOAT(IOOT) 

2 CONTINUE 
LANE=255 
HALFLN=128, 

sa«=o 

DO 6 1=1,100 
SUK=SUK + X (I) 

6 CONTINUE 
AVE=SUN/100. 

DO 7 1=1,100 
X(I) = X{I)-AVE 

TF(X(X) .GT.HALFLN) X (I) =X (I) -LANE 
IF(X (I) .LT.-HALFLN) X (I) =X{I) +LANE 
C CHANGE UNITS FROM 255THS TO CEC. 

X(I) =X(I) *100,/ 

7 CONTINUE 
DO 3 1=1,900 
CALL DATA(IOUT, 

3 CONTINUE 
DO 4 1=101,200 
CALL DATA(I0UT, 

X (I) = FLOAT (lOUT 

4 CONTINUE 

sun=o. 

DO 8 1=1,100 

suK-sua+x{ioo*i 

8 CONTINUE 
AVE=SUM/100. 

DO 9 1=101,200 
X(I)-X{I)-AVE 
IF(X (i> .GT.HALFLN) X{I)=X(I) -LANE 
IF(X(I) .LT.-UALFLK) X (1) =X (I) *LANE 

C CHANGE UNITS FROM 255THS TO CEC; 

X(I)=X(I) *100./255. 

9 CONTINUE 

C CALL AUTCOP. SUBROUTINE TO COMPUTE AND PLOT R (T) . 

CALL AUTCOR(X) 

PRINT 103 

103 FORMAT (lx, ' THIS IS THE AnTOCORBELATlON L'F REAL OMEGA NOISE*) 
PRINT 104, I SKIP 

104 FORMAT (IX, 'THE DATA STARTED ', 15 ,» SAMPLES AFTER THE STABT OF 
CTHE D TIME SLOT') 

C GALL THE HISTOGRAK SOF-EOUTINE FOP REAL OMEGA NOISE PHASE . 

CALL HIS TO (X) 

PRINT 106 

106 FORMAT (IX, 'THIS IS THE HISTOGRAM OF REAL OMEGA NOISE PHASE VS. 
C PHASE IN CEC.') 

PRINT 104,ISKIP 
100 STOP 
END 


255. 

S1G0) 

6100) 


) 


) 




Figure E-13. (Continued) 
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xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx. 

SObROUTlNE AUTCOS(X) 

EEAL X(200) ,XM 100 ),M 100 ),XXA( 100 ),XXL 12) ,vn(12) 

DATA YYL/*AS*tl», t, r I' » ^ ^ / 

COaPUTE AND PLOT AUTOCORRELATION OF THE NOISE, 

DO 10 1=1, 100 
XX (I) =X(I) 

10 CONTINUE 

; CORPUTE AUTOCOBFELATION FOR ONE VALUE OF TAO. 

DO 20 11=1,100 
SUM=0. 

DO 30 1=1,100 
IF(I1.20. 1) GO TO 16 
GO TO 15 

16 J=I 

GO TO 17 
15 J=I+II 

17 CONTINUE 
PROD=XX (I) «X (J) 

SUN=S0a+PROD 

30 CONTINUE 

R(II) =SUPI/100. 

20 CONTINUE 

DO 50 1=1, 100 
XXA (I)=FLOAT (I) 

50 CONTINUE 

CALL NKPLOT (XXA,R, TOO, 1 ,XXL,YYL) 

PETUEN 

G XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SOBKOUTINE GR A NGE ( Y, INT VL , A ,SH ALLA , TNF.) 

INTEGER TME 

REAL Y(500 0), INTVL 

TA0=INTVL 

E=EXP (^SKALLA*TA'J) 

PSI= (A*A*SMALLA) 
coKs=sQ8T(Psi* (1.- (R*n) ) i 
CALL ANOBM(AN) 

Y{1)=AN*PSI 
DO 10 1=2, TNE 
J=T-1 

CALL ANORH(AN) 

9 FOnflAT(lX, 1OE10. 2) 

Y(I) =AN*GON3<-R*y (J) 

10 CONTINUE 
return 

e SxxxxxxxxxxHxxxxxxx*xxxxxxxxxxvxxxxxxxxxxxxxxxxxxixxxxxxxxxxxxx- 

ONE OATA SABPI.E TO »At» MOGAA, 


G..,ON TAPE EOF, RETURNS TO « 

c!!* R. W. LILLEY, AVIONICS, JANUARY, 1976 

Figure Et- 13, (Gontinued) 






n n n n o r. n rj r. o rj n 


INTEGER 1/0/ 

EOGICAL^l IN (1000) 

EQUIVALENCE (II, INE(1) ) 

11 = 0 

1F(I. EQ.0)GO TO 10 
30 INB(4)=IN(I) 

I0UT=II 

1=1*1 

IF (I. GT, 1000) 1=0 
BETUBN 

C • • • 

10 1 = 1 

17 READ{11, 5, END = 2G) IN 
5 FORMAT (5 (200A 1) ) 

GO TO 30 

C» • • 

20 RETURN 1 

C. .. 

ENTRY START 

...STARTS DATA OUTPUT AT A RANDOM POINT ON TAPE. READS OP TO 
25 RECORDS, AND STARTS BETWEEN BYTE 1-1000. 

...SEED KITH UNDEFINED VAxRIABLE 
IF (IQ. Eg. 0) IQ=5**1 3 
IQ=IADS (IQ/2) *2*1 
CALL INTGEN(IQ) 

RETURN 

ENTRY RAND{NBYTE) 

GALL IPAND(NREC, 24) 

GALL IKAND(NBYTS,999) 

... 

DO 51 IR = 1,NREC 

51 READ (11, 5) IN 
I=NBYtt* 1 

PRINT 52,NKEC,NBYTE 

52 FORMAT (• THIS RUN STARTS AFTER */lX,I2,* SKIPPED RECORDS AND*, 
♦1X,I4,‘ SKIPPED BYTES •) 

GO To 17 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE ILPF (Y , T AU , A , 0M , NUM) 

PRINT 14,TAU,DW 
14 FORMAT(1X,2F10.3) 

EOUTiNL RETURNS »NUK* IDEAL LOW-PASS-FILTERED, GAUSSIANLY 
DISTRIBUTED RANDOM VARIABLES. 

Y IS THE VECTOR OF SIN(X)/X GORRELATED NUMBERS; TAll IS THE INDEP- 
ENDENT TIME SHIFT VARIABLE; A IS A HASTE; BW IS THE BANDWIDTH 
OF THE FILTER (ONE-SIDED) ; NOH IS THE NUMBER OF NUMBERS HETURNED. 
REAL Y(5000) 

H= {S I N ( 6 , 2 B 3 1 B5* BK *T A U) ) / (S . 28 3 1 8 5*= BW *T A U ) 

CONS = SQRf (ABS (1.- (R*R)) ) • 

PRINT 12, R, CONS 
12 FORMAT (1X,2F1 0.3) 

GALL ANORn(AN} 

Figure E-13. (Continued) 
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y ( 1 ) =AN 

DO 10 1=2, NUM 
J=I-1 

CALL ANOI?M(AN) 

Y(I) =AN«CONS + 8*y (J) 

11 F0RHAT(1X,2F8.2) 

1G CONTINUE 
KETURN 
BSD 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE HISTO(X) 

G ROUTINE PLOTS HlSTOOPA« OF NUKBERS IN VECTOR (X(200)). 

INTEGER IX(20G) 

HEAL X(200), XA(1G1), YA(101), XL(12), YL(12) 

DATA XL/*N» ,«0' ,*I» ,‘S* ,» E* ,* * , * P* , * H' , » A • , * S • , • E* , • •/ 

DATA YL/'H*,'I*,'S* ,*T*,'0* ,*R','R»,»A*,*K',' • • •/ 

DO 10 1=1,200 
IX(I)=INT(X(D) 

10 CONTINUE 

C CREATE X-AXIS VECTOR OF -50. TO 50, CEG. 

DEV=-50. 

DO 60 1=1, 101 
YA (I) =0. 

XA (1) =DEV 
DEV=DEV+1. 

60 CONTINUE 

C CREATE Y-AXIS VECTOR Of FREQ. OF OCCURENCE OF FACE CEC VALUE, 

DO 70 1=1,200 
DO 80 vl= 1,101 

IF(IX (I) .EQ.INT(XA (J) ) ) YA (J) =YA (J) f 1, 

00 CONTINUE 
70 CONTINUE 

C NOEHALIZE RELATIVE FREQ. OF OCCURSNCE. 

DO 90 1=1,101 
YA(I)=YA{1) /200. 

90 CONTINUE 

CALL NMPLOT (XA,YA, 101,1 , XL, YL) 

HETOHN 

END 


Figure E-1 3, 


(ConMnued) 
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OMEGA SENSOR 
MAIN PROGRAM 


PROCESSOR SIMULATION, 


G 

C 


G 

C 


G 

C 

G 


F. SALTER, AVIONICS ENGINEERING CENTER, MARCH 1976 

«4i ♦ ♦ ♦ ♦ ♦ ♦ * * * ♦ • ♦ • • * • • ****♦♦< 

OMEGA AUTO-SYNC ROUTINE 
integer lop (1 GOG) 

INTEGER X(1G00) , DIEF (1C00) , C,LI VSLT ( 100G) ,CORR,CTR ,SYNCPT 
INTEGER SLOT (1000) 

INTEGER SMDIFF (1G00) 

INTEGER STA{8) 

INTEGER WITHIN, V (1000) 

INTEGFH IMIIR(IOOO) 

INTEGER PRASE (3) 

INTEGER XX (2043) ,YY (2048) ,XL(12) ,YL(12) 

INTEGER* 2 TSLOT (8) /9 0, 1 00, 1 10, 1 2 0 , 90 , 1 1 0, 1 20, 1 00/ 

INTEGER LCW(8),PHI 

INTEGER CPHI(8,10G0) ,CCPHI(1000) 

INTEGER II (8) ,SOGAIN (8) ,VEL (8) ,RATB 

INTEGER CONST, 3UMDIF (8) ,V AR (8,1000) , WAR (1000) 

DERIVE random’ STARTING POINT ON DATA-BASE TAPE 
SKIP •NSKIP* FILES ON DATA BASE TAPE. 

NSKIP^U 

CALL START 

CALL RAND(NBYTE) 

TAKE GNE RECORD (1000 BYTES) OF DATA AND HIRROP II 

DO 1 1=1,1000 

CALL DATA(X (I) ,ft50) 

CALL WIRR (X (I) ,INMIR (I) ) 

1 CONTINUE 

CALL PL0T1G(X,1) 

CALL PLOT10 (INMIR, 1) 

OBTAIN A 1000 BYTE VECTOR, OF DIFFERENCES TO BE USED AS A MEASURE 
OF PHASE COHERENGP IN EACH 10 MSEC SAMPLE SLOT 
CALL DIFFER (INMIR, DIFF) 

CALL PLOT 10 (DIFF, 1) 

4 F0RMAT(1X,2SI4) 

SMOOTH DIFFERENCES (OPTIONAL), 

CALL SMOOTH (DIFF, SflDIFF) 

PRINT 151 

151 F0RKAT(1X,»9 PT. STRAIGHT AVERAGING SMOOTHER INSTALLED*) 

GENERATE A 1000 BYTE VECTOR OF V'S AND 0*S REPRESENTING SIGNAL 
AND NOISE SLOTS 

MUST SET THRESHOLD "C” TO COMPARE DIFFERENCES TOi 
G= 1 20 

12 CONTINUE 

CALL eOMPAP (G, SMDIFF, LIVSLT) 

PERFORM BILEVEL CORRELATION 
eORR=0 
J=0 

CALL PATEEN(SLGT) 

10 CALL BICOR (LIVSLT, CTE,SLOr) 

1= J ♦ 1 

IF(GTB.GE,CORR) GOER=CTR 


OF KNOWN OMEGA PATTERN TO LIVE DATA 


Oi? 




E^14. FORTRAN Program RSOSP Used to Simulate Autpmat'c Synchronization 
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IF(CTK.£Q.COr,R) SYNCPT=J 
J=J*1 

1F(J.GT. 1000) GO TO 100 
NSHII T=1 

CALL 3MIFT (LIVSLT,NSHIFT) 

GO TO 10 

C PRINT CORRELATION VALUE AND SYNCPOINT VALUE OBTAINED. 

100 PRINT 110,C,COPR,SyNCPT 

110 F0SNAT(1X,'C = ',I3»3X,« CORR = * , 14 , 3X, ‘SY NCPT = *,I4) 

«ITHIH=10* ( 330+NDYTE-SYNCPT) 

PRINT 13, WITH IN 

13 FORMATdX, 'SYNC W ITHIN* , 15 , 1 X , » NSHC • ) 

NSHIFT^SYNCPT 

C SHIFT VECTOR 0? PHASE VALUES TO SYNCPOINT AND PLOT THEM OUT. 

CALL SaiFT (X, NSHIFT) 

CALL PLOT10 (X , 1) 

GO TO L>0 
PRINT 111 

111 FORMAT(1X,'ALL STATIONS ON AIR’) 

C RESET DATA BASE TAPE TO SYNCPOINT. 

NSniFT = 1 000-NSllIPT 
DO 14 1=1, NSHIFT 
CALL DATA (IOUT,i>50) 

14 CONTINUE 

C OMEGA BURST FILTER ROUTINE. 

C PEEFOEM BUF.ST PLL FILTERING TO INCOMING OMEGA PHASE MS ASU RFM ENTS , 

DO 18 1=1,8 
LCK (I) =0 
SOGATN (I) =0 
VEL(I)=0 

18 IX (I) =1 

C READ 1000 OMEGA FRAMES. 

DO 26 M= 1,1 000 

C READ a TIME SLOTS IN EACH OMEGA FRAME. 

19 DO 25 1=1, a 
JJ = TSLOT(I) 

C READ TSLOT(I) PHASE MEASUREMENTS IN EACH TIME SLOT. 

LPGTR=C) 

PATE=0 
SUMDIF (I) =0 
DO 20 J=1,JJ 
CALL DATA (PHI, S50) 

CONST -2 

CALL 5MAPLL(PHi,LPCTS,LCN(I) , CONST, RATS, SUMDIF (I) ) 

IDIF= SUMDIF (1) 

LLCW=LCW(I) 

21 FORMAT (IX, 5110) 

20 CONTINUE 

22 FOEMAT(IX) 

CPHI(I,iI (I) ) =LCK (I) 

C COMPUTE SIGNAL VAFIANCE IN TIMESLOT. 

SUHDiF(I) =SIJMDIF (I)/TSLOT (I) 

VAR|1,II (I) ) =SUM0IF(I) 

IF (HATE.LT.O) VEL{I) =VEL{1) “EATE*SOGAIN(I) 


Figure E-14. (ConHnued) 


c 


IF(RATK.f;T. 0) VEL (I) =VF.L (I) + E AT E*SOr, ftl N (I) 

Lew (1)=LCW (I) ♦VEL{I) 

IF(LCW(1) .LT.O) LCH (I) =256+LCH (I) 

IF (Lew (I) . GT. 255) LCK (I) =LCM (I) -256 
II (1) =il (I) >1 

SKIP SAMPLED FROM ,2 SE'C. NOISE GAP. 

DO 2J K=1,20 
23 CALL DATA (PHI, «50) 

25 CONTINUE 

26 CONTINUE 

C UUCP EACH STATION’S TIHESLOT PJIASS ESTIMATE INTO VECTOP TO BE PLO- 

C TTED BY PLOT1 SUBROUTINE. 

30 DO 40 1=1, B 
DO 3*^ J = 1,10f)0 
CCPlil (0) =CPHI (I, J) 

C DUMP EACH STATIONS VAPIANCE INTO A VECTDP TO BE PRINTED OUT. 

WAR (J) =VAP (I, J) 

35 CONTINUE 

C CM EG A LOP ROUTINE 

PRINT 3e),l 

36 FOFMATC PHASE OF STATION ’,12,* MINUS CLOCK') 

PRINT bO,GCPHl 

60 FOFMAT (1 X, 1 0110) 

PRINT 61,1 

61 fORMATC SIGNAL VARIANCE OF STATION ',12) 

PRINT 6 0 , WAR 
GALL PLOT1 (CCPUI, 1) 

40 CONTINUE 
IS1A1 =3 
ISIA2=« 

CALL LOrO (CPU 1 , 1ST A 1 , IS TA2, LOP) 

G PLOT LOP SPECIFIED BY ISTA1-ISTA2. 

CALL PLOT1 (LOP, 1) 

STOP 

50 STOP 50 
END 

SUBROUTINE DIFFEE': (X , DIFF) 

C OMEGA auto-sync SIMULATION 

C 

C ROUTT SL TAKES 1000 BYTE VECTOR OF INPUT PHASE DATA AND RETURN 

C DIFFERENCE VALUE WHICH IS INVERSELY PROPORTIONAL TO THE PHAS 

C COHEKENCE (SIGNAL STEENGTH) IN EACH 10 MSEC SAMPLE SLOT. 

C 

C R. SALTER, AVIONICS ENGINEERING CENTER, MARCH 1976 

C 

INTEGER DIFF (1000) , X (1000) 

DO 20 1=1,9 99 
J=I + 1 

DIFF (I) =IABS (X (I)-X(J) ) 

20 CONTINUE 

DIFF (1000) =IABS (X(1) -X (1000) ) 

RETURN 
END 

SUEROUTISE GOMPAE (G , DIFF, LI VSLT) 



(GonHnued)' 


Figure E-14. 
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C OMEGA AUTO-SYNC SIMULATION 

C 

C ROUTINE COMPARES 1000 BYTE VECTOR OF DIFFERENCES TO THRESHOLD 

C CONSTANT AND RETURNS A BILEVEL VECTOR CONTAINING t'S FOR SIGNAL 

C AND 0*S FOR NOISE 

C 

C R. SAi.TLR, AVIONICS ENGINEERING CENTER, MARCH 1976 

C 

INTEGER C,DIFF(1OO0) ,LIVSLT (1000) 

DO 10 1=1,1000 
IF(DTFF(I) .GE.G) LIVSLT(I) =0 
IF (DIFF(I) .LT.C) LIVSLT (I) =1 
10 CONTINUE 
RETURN 
END 

SUBROUTINE BI COR (LI VSLT ,CTR , SLOT) 

C OMFGA AUTO- SYNC SIMULATION 

C 

G ROUTINE PERFORMS BILEVEL G HOSS-CORR ELiri ON OF OMEGA TRANSMISSION 

G PATTERN OFF AIR WIT!! THE KNOWN PATTERN GENSRATED BY "PATTERN** , 

C AND RETURNS VALUE PROPORTIONAL TO CORRELATION 

G 

C lu SALTER, AVIONICS ENGINEERING CENTER, MARCH 1976 

C 

INTEGER LIVSLT (1000) ,GTR 
INTEGER SLOT (1000) 

CTF. = 0 

DO 10 1=1, 1000 

IF (LIVSLT (I) .EQ.SLOT (I) ) CTH-GTR + 1 
10 CONTINUE 
RETURN 
END 

SUBROUTINE SHIFT (LIVSLT , NSHIFT) 

C OMEGA AUTO-SYNC SIMULATION - 

C 

G ROUTINE SHIFTS LIVSLT VECTOR "NSHIFT" SLOTS TO THE RIGHT. 

C 

C R. SALTER, AVIONICS ENGINEERING CENTER, MARCH 1976 

C MODIFIED BY B. SALTER, APRIL, 1976 

C 

INTEGER LIVSLT(IGOO) ,AUK5LT (1000) 

DO JO K=1, NSHIFT 
DO 10 1=2,1000 
J=l-1 

AUXSLT (I) =LIVSLT (J) 

10 CONTINUE 

AUXSLT (1 ) =LIVSLT( 1000) 

DO 20 1=1,1000 
LIVSLT (I) =ArJ':SLT (I) 

20 CONTINUE 
30 CONTINUE 
PETURH 
END 

SUBROUTINE PATSEN(SLOT) 

C OMEGA AUTO-SYNC SIMULATION 


Figure E-14. (ConHnued) 
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c 

c 

G 
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ROUTINE CREATES KNOWN OMEGA TRANSMISSION PATTERN 
B. SALTER, AVIONICS ENGINEERING CENTER, HABCH 1976 
INTEGER SLOT (1000) ,STA (8) 

INTEGER Jj <6) /I, 111 ,231 ,361 ,501,631,741,88 1/ 
INTEGER KK (8) /90, 2 1 0, 34 0, 48 0, 6 1 G ,720 , 860 ,980/ 
READ(1,1)STA 
1 FOEMAT(8I1) 

DO 10 1=1,8 
J=JJ(I) 

K=KK(I) 

rOO DO 110 L=J,K 
SLOT CD =STA (I) 

110 CONTINUE 
H=K+1 
N=K + 2 0 

DO 12 0 L=K,S 
SLOT (I.) =0 
120 CONTINUE 
10 CONTINUE 
130 RETURN 


END 

PSTURSiS ORE DATA SAdPlB TO BAIS PROGRAl* 


!!!oN tape EOF, RETURNS TO * 


P. W. LILLEY, AVIONICS, JANUARY, 1976 
INTEGER 1/0/ 

LOGICAL^I IN (1000) ,INB (4) 

EQUIV ALENCE(II,INB(1) ) 

11=0 

IF (I.EQ.O)GO TO 10 
30 INfi(4)=IN(T) 

IOUT=II 
1=1 + 1 

IF {I.GT. 1000) 1=0 
EETUEN 

C ■ ■ « 

10 1=1 

17 READC11, 5,tND=20) IN 
5 FORMAT (5 (2C0A1) ) 

GO TO 30 

G« • • 

20 RETURN 1 . 

c • • • 

ENTRY START 

e.i. STARTS DATA OUTPUT AT A BANBOH POINT ON TAPE. 
G... 25 RECORDS, AND STARTS BErWEEN BYTE 1^1000. 

c» • • 

g!!*.SEED with UNDEFINED VARIABLE 
iF(IO.EQ.O) IQ=5»*13 
IQ=IABS (IQ/2) *2+1 


READS OP TO 


Figure E-14. (GonHnued) 
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CALI I.NTGI:N(IQ) 

PBTUfiN 

ENTRY HAND(NnYTE) 

CALL TPAND (MREC,24) 

CALL ir<AND(NBYTE,999) 

C • • • 

DO 51 IR = 1,NBEC 

51 READ(11,5) IN 
I=KbYTE>1 

paiNT 52, NFEG, MBYTE 

52 FOBHATC THia RUN STARTS AFTER SKIPPED RECORDS AND', 

♦1X,14,* SKIPPED BYTES *) 

GO TO 17 
END 

SUBFOIJTINE MIRR(IK,INMIH) 

C. ..OMEGA DATA SIMULATION 

C... ROUTINE RLFLKCTS COUNTS 123-255 MIFRORED TO 127-0 

C.y, R. K. LILLEY, AVIONICS, JANUAPY, 1376 
C nODlFILD BY P. SALTER, APRIL, 1976 

IF(IH.LE.127) GO TO 10 
INMIh=255-IN 
RETURN 
10 INM1E=IK 
REIURN 
END 

SOBRCUTINE SMOOTH (DI E F, SMDI FF) 

C OMEGA AUTO-SYNC SIMULATION 

C 

C ROUTINE SMOOTHS VECTOR OF 1000 DIFFERENCES DY PASSING EACH SAMPLE 

C POINT THRU A 9 POINT EXPONENTIAL FILTER 

C 

C F. SALTER, AVIONICS ENGINFERING CENTER, MARCH 1976 

G 

INrEGSR SMDIFF(IOOO) 

INTEGER DIFF(IOOG) 

DO 10 1=1,1000 
IF (I. LT. 5) GO TO 20 
TF(I.GT.095) GO TO 30 
J=I-1 
K=T-2 
L=I-3 
M=I-4 
25 JJ=I*1 
KK=I+1 
LL-H^3 
MH=I+4 

C TEMPORARY SMOOTHER NINE PT. STRAIGHT AVERAGER. 

C OSES 4 RTS. BEHIND , 4 ?TS, AHEAD, AND CURRENT PT. EQUALLY WEIGHT 

35 SMDIFF (I) -DIFF(I) +DIFF(J) +DIFF (K) +DIFF (L) ♦DIF? (M) ♦DIFF(JJ) ♦DIFF 
♦ (KK) +DIFF (LL) +0IFF (MM) 

10 CONTINUE 
GO TO 40 

20 j= 100 0- MOD ( 1 00 0-^T^ 1 i 1 00 0) 

K=10QO-HOD (1000^1+2,1000) 


Figure £•^14, (Continued) 






o Q n o 


L=1000-KOD(1000“I*3, 1000) 

K=1000-HOD (1 000-14- U, 1000) 

GO TO 25 
30 J=l-1 
K*I-2 
L=I-3 
W=I-4 

JJ=MOD (14-1, 1000) 

KK = MOD(l4-2, 1000) 

LL=riOD (14-3, 1000) 

KH=ROD (14-4, 1000) 

60 TO 35 
40 BETUHN 
END 

SUBSOOTINE PLOT10 (5C,LLL) 

C BOUTINE PLOTS DATA VECTOE ”X” 10 POINTS PER LINE HITS STABS.,, 

C 

C SCALE OF PLOT IS 1-128,, . IP VALHSS OF "X'VARE LESS THAN 128 THEN 

G OPERATOR SHOULD SET LLL=0 AND A STRAIGHT PLOT HILL RESULT, IF "X" 

G HAS VALUES BETWEEN 128-256 LLL SHOULD BE SET TO 1 AND SOUTINE WILL 

C HALVE EACH VALUE BEFOKE PLOTTING, 

C 

C R. SALTER, AVIONICS ENGINEERING CENTER, APRIL, 1976, 

C 

INTEGEE^ BK/‘ '/, ST/* ♦ V, LIN (128) , X ( 1 000) 

K=0 

DO 30 J= 1,991 
DO 10 1=1,128 

10 LIN(I)=BK 
J=K+ 1 
K=J4-9 

DO 20 I=J,K 
IF (LLL.EQ.O)GO TO 12 
KN=X(1)/24-1 
GO TO 13 

12 CONTINUE 
«M=X (I) 

13 CONTINUE 
LIN(MN) =ST 

20 CONTINUE 

PRINT n,LIN 

11 FOP»1AT(1X,128a1) 

30 CONTINUE 

RETURN 
END 

SUBROUTINE LOPO (CPHI , ISTA 1 , IST A2, LOP) 

ROUTT NECSEATES A VECTOR OP LOP POINTS (LOP(IOOO)) BY DIFFSEENCING 
STA1^3TA2 PHASE REASUPEMENTS SUPPLIED IN GPHI(8,10OO) ARBAY. 

B. SALTER, AVIONICS ENGINEERING CENTER, OHIO UNIVERSITY, JULY 1976 

INTEGER LOP (1000), CPHI (8, 1000) 

DO 10 1=1,1000 

LOP{I) =CPHI(ISTA1,I) -CPHI (ISTA2,T) 

IF(LOP(I) .LT.0) LOP(T) -LOP(I) 4-256 
10 CONTINUE 


Figure E-14. (Gontinued) 


no n o o n n n n ft o 


fiETURN 

END 

SOEROIJTINE SMAPM. (PHI, LPCTH,PHIREF, CONST, BATE, SUNDIF) 

ROUTINB PEKFORHS DIGITAL FItTEBING GN INCOflING OMEGA PHASE 
MEASUREMENTS ( PLL, BUEST FILTER). OPERATION AND riHE CONSTANTS 
CLOSELY RESEMBLE HARDWARE MAPLL. 

* ’CONST** IS A VALUE TO BE SPECIFIED FOB BAX. LOOP CTB. VALUE 
BEFORE ADJUSTING PHIKEF. (LOOP TIME C0?1STANT OF AMOUNT OP INTEGRA. 

B. SALTER, AVIONICS ENGINEERING CENTER, JULY 1976. 


C 

C 

G 

C 

G 

G 


INTEGER PHI, PHIREF, CONST, B ATE, LPCTR 
INTEGER DIFF,SUMDIF 

COMPARE INCOMING PHASE TO LOOP CONTROL WORD (LCW) (PHASE REFERENCE 
DIFP^PllI^PniRSF 

KEEP A RUNNING SUN OF ABSOLUTE VALUES OF DIFFERENCES TO INDIGATP 
SIGNAL VARIANCE, 

IF(DIEF) 10,50,20 

CHECK FOR EDGE OF LANE AMBIGUITY, 

10 IF (DIF?, LT, -128) GO TO 25 
SUMDlF=SUMDlF + rABS (DIF?) 

GO TO 40 

20 IF (DIFF.GT. 12 8) GO TO 35 
SUMDIF=SUMDIF+rA3S (DIFF) 

GO TO j6 

INCREMENT LOOP COUNTER. 

30 LPCTR=LPCTR+1 

EXAMINE LPCTR FOR INCE, MORE THAN ’CONST* TIMES CORRELATED. 

IF (LPCTR .LE. CON ST) RETURN 
LPCTR=0 

PHIREF-PHlP.EFtl 

HATE=RATE+1 

CHECK FOR MODULO 255 LANE COUNTS, 

IF(PHIPEF.GT. 255) PHIREF=0 
RE TURN 

DECREMENT LOOP COUNTER. 

40 LPeTP=LPCTE- 1 

EXAMINE LPCTR FOR DECS. MORE THAN 'CONST* TIMES CORRELATED, 

IF (LPCTR. GE. -CONST) RETURN 
LPCTE=0 

PHIREf=PHIREF-1 


RATE=SATE-1 

CHECK FOR MODULO 255 LANE COUNTS. 
IF (PHTREF.lt. 0) PHIEEF=255 
50 RETURN 
25 PHl=PHI*-256 

DIFF=PHI-PHIRSF 
SHMDIF=SUMDIF+IABS (DIFF) 

GO TO 30 
35 PHI=PH1^256 

DIFF=PH1-PKIHEF 
SUMDIF=SUMDIF+IABS (DIFF) 

GO TO 40 



END 

SUBROUTINE PL0T1(X,LL) 



EOOTISE PLOTS DATA VECTOR '‘X" WITH STARS, ONE POINT PER LISE ON A 
1-12B SCALE. 


C 

c 
c 

C IF VALUES OP "X" CAN FALL BETWEEN 128-255, SET LL=1 AND ROUTINE 

C WILL HALVE EACH VALUE BEFORE PLOTTING. OTHERWISE, IF VALUES OF 

C "X« ARE ALL LESS THAN 128, ROUTINE KILL PLOT STRAIGHT VALUES IF 

C LL-0 

C 

C MUST DIMENSION X = VECTOR OF VALUES TO BE PLOTTED, 

C E.W. LILLEY, AVIONICS, JANUARY, 1976 

C MODIFIED BY R. SALTER, AVIONICS, APBrL,1976. 

C 

INTEGER BK/’ •/» ST/* ♦ V,LiN (128) 

INTEGER X(1000) 

DO 10 1=1,128 

10 LIN (I) =BK 

DO 20 1=1,200 
IF(LL.BQ.0)GO TO 12 
MM=XiI) /2^1 
GO to 13 

12 CONTINUE 
KM=X(I) 

13 CONTINUE 
LIN (14 M) =ST 
PRINT 11, LIN 

11 FORMAT( IX, 128A1) 

LIN(MM)=BK 

20 CONTINUE 
EETUI'N 
END 


Figure (ConHnged) 
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SMAPLL SIKULATIOK ANALYSIS TO OPTIMIZE LOOP PABAMETEBS 
B. SALTER ; JANUAFY 1977 ; AVIONICS ENGINEERING CENTER. 

REAL GA{5000) ,GB{5G0G) ,G(50G0) ,INTVL,XX (100) ,YY (10G) 

BEAL S {f)00) ,STD (600) 

INTEGEP NN (5000) ,ESrtOR (200) 

INTEGER NOISE (5000) ,TME,PHI (5000) ,PH1IN, PHIOUT ,SUMDIF,CTR2 
INTEGEh E3ItSUM,TSLOT, PHASE, DELTA , CONST, RAT R,HALFLN 
DIMENSION YYY (101) ,XXX(K'1) ,XXXX (101) ,YYYY (101) 

BEAL XXXL(12) ,YYYL (12) ,XXXXL(12) , Y YYYL ( 1 2) , XXL ( 1 2) ,YYL(12) 

DATA XXL/'S* , ’I^'G* A*,»L» , V'f 

DATA YYL/»P»,»L',<L',* */ 

XXXXXXXXXXX •'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


INITIALIZE PARAMETEHS: A=5IGNAL APLITOOE; LEVELS =NUSBER 
LOOP QUANTIZING LEVELS; CONST-NUMBER BITS SKIPPPED IN LOOP 
BIDIRECTIONAL COUNTER (LOOP FILTER); INTVL=SAKPL ING RATE IN 
SECONDS; VEL=VFLOCITY OP PHASE RASP INPUT TO LOOP; PHASE- 
INITIAL PHASE INPUT; SOGAIN=SECOND ORDER GAIN. 

CONTRL IS A CONTROL WORD TO SPECIFY OPTIONAL GA'LCULATIGNS AND PLOT 
CONTEL=0. 

LL=1 

A=3. 

S0GA1N=. 00 
VEL1=300. 

VEL^VELI 

PRINT 113,SOGAIN,VSL 

113 FQRMAT(1X, "SECOND ORDER GAIN = * ,F5 . 2 , 1 OX , • PH AS E VELOCITY = * , F7. 2) 
HALFLN=64 
LANE= 127 
LEVEL5=120 
TSLOT=1UO 
MUS=5000/TSLOT 
CONST-1 
SMALL A=1B9. 

INTVL=^.01 

TME=5000 

C A P A = 2 . * SQE T ( 4 , /S M A LLA) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

GENERATE NOISE VIA CALLS TO RANDOM NUMBER GENERATORS AND 
LOW PASS FILTERS TO CORRELATE IT, 

CA LL I NT GEN (5**15) 

CALL GRANGE (G A, I NTVL, CA P A,S M ALL A,T ME) 

COMPUTE THE MEAN AND VARIANCE OF NUMBERS PETURNSD FROM RANDNR, 

AND THROUGH RC LOW PASS FILTER GRANGE. 

CALL MNVAR(GA, AVE, VARI) 

CALL BNTGEN (5**1 1) 

CALL GnANGD(GB,INTVL,CAPA,SMALLA,TME) 

9 FORMAT (IX, 1GE10.2) 

SNK=10.*ALOG1Q(A*A/(VAFI*2.) ) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


Figure E^IS. FORTRAN Program RSSIM Used to Simulate SAAAPLL2 and Evaluate 


Performance. 


-V08- 


nn o onn nnn 


PEKFORK CALCULATION OF LOOP VARIANCE FOR SEVERAL SNRS 
SEVERAL TIMES. EACH TIMS INCREASING THE LOOP TIME CONSTANT 
(M) (MORE INTEGRATION OR SKIIPED BUTS IN THE COUNTER). 

DO 45 MM=1,8 
PRINT 46,SOGAIN 

46 F0EMAT(1 X, 'SECOND ORDER SAIN = »,F5,2) 

DO 40 H=l, 1 
CONSI = T6 
PRINT 11, CONST 

11 F0BMAT(3X, • LOOP CONSTANT = *,I5) 

PERFORM CALLS TO SMAPLL AND CALCULATE OUTPUT VARIANCE FOR 
100 GIVEN SNRS , INCREASE SNR EACH AROUND THIS LOOP. 

DO 30 L=1, 1 

CALCULATE KE« NOISE PHASE INPUTS FOR EACH 3MP-. 

3 DO 1G 1=1,5000 

G(I) =ATAN2 {GR (I) , (GA (1) *A) ) 

G(I) =G (I) *57. 2057795 
114 FORMAT (1X,E9. 1) 

112 FORMAT (1X,E12.3) 

NOiSS(I) =INT{G(1) ) 

10 CONTINUE 

IF (CONTRL.LT. 10.) GO TO 14 

CALL lilSTi subroutine TO PLOT HISTOGRAM OF INPUT NOISE PHASE. 

CALL iilSTi (NOISE, 5000, SNR, VAR) 

14 CONTINUE 
K=TSLOT 
ERPSUM-0 
SUMDIF=0 
GTR2=0 

calculate rate for 2ND ORDER LOOP STEADY STATE TRACKING FOP GIVEN 
VELOCITY SIGNAL INPUT. (I.E. START LOOP OUT BEING LOCKED IN SS) . 
IF(SOGAIN.EQ.0.) GO TO 17 

RArE=INT(VLL*FLOAT(L£VFLS) /(5760,*SOGAIN) ) 

17 CONTINUE 
PHIOUT-0 
PAT‘F=0 
LPCTR=0 
JJ=1 
J=0 

VEL=VEL1 

DELTA=INT (VEL/16.) 

PHAS E=0 

C CALCULATE PHASE INPUT TO LOOP BY ADDING TRUE PHASE PLUS NOISE PHAS 

15 DO 20 1=1, K 
PHI (I+J) -PHASE 

* PHIIN=PHI(I'-J) +NOISE (I+J) 

PEiriN=MOD (PHIIN, 360) 

IF{PHIIN.LT.O) PHIIN=PHIIN + 360 
] X=FLOAT (PHilN) 

CALL QUANT(X, LEVELS) 

PHriN=INT(X) 

PIIIIN=MOD (PHIIN, LEVELS) 

IF(PHIIN.LT.0j PHIIN=PHiIN*LANE+1 

CALL SKAPLL(PHIlN,LPGTE,PHIOUr,CONST,CTR2,3nHDIF) 
in FORMAT (IX, 614) 


Figure E-15. (CbnHhued) 
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20 GOSTIKUE 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKiCXXXXXXXXXXXXXXXXXXXXXXXXXXX 

CALCtlLATS LOOP ERROR BY TAKING PHI (TRUE PHASE IN) - PHIOUT 
(LOOP estimate of PHI) . 

X=FLOAT (PHI (I+J) ) 

CALL yUANT (X, LEVELS) 

PHI(1+J) =INT (X) 

PHI (I+J) =H0D (PHI (I + J) , LEVELS) 

IF(PHI(I+Jj .LT,0)PHI(I + J) =PHI (I+J) +LANB+1 
IQUANT=PHI (I + J) 

LQUAKT=PHIOUT 
C0NTRL=1 1. 

IF (CONTRL.lt, 10.) GO TO 1ft 
CALL PLTPHS (LQaANT,IQUANT) 

16 CONTINUE 
CONT?L=0. 

ERROR (JJ) =PriIOUT-PHI (I+J) 

IF (ERHOll (JJ) . GT, HALFLN) ERROR (JJ) =£ROOR (JJ) -LEVELS 
IF (ERROR (JJ) . LT.-HALPLN) EPROR (JJ) =RRROS ( JJ) +LEVELS 
ERROR (JJ) =INT (FLOAT (ERROR (JJ) ) +5.625) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxyxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PERFORM RATS^ AIDING FOP 2ND ORDER SHAPLL, 

RATE=RAT E+CTR2 

XA = FLOAI (PHIOUT) +FLOAT(EATE) +S0GAIH 
PH10UT=INT (XA) 

PHIOUT=HOD ( PHIOUT, LEVELS) 

IF (PHIOUT. LT.O) ?HIOUT=rilIOUT + LEVLLS 

CALCULATE PHASE AT TIME OF NEXT TIME SLOT (10 SEC LATER). 
DEI.TA=INT (VEL/16.) 

PHASE=PH ASE+DELTA 
PHAS E=MOD (PHASE, 36 0) 

IF (PHASE. LT.O) PH ASE=PHA5 F+ 3 60 

J=J+TSLOT 

SUMDIF=0 

CTR2=0 

LPCTF-0 

JJ=JJ+1 

I?(J.r,T.2500) VEL=-^VEL1 
IFlJ.LT,5000) GO TO 1 5 

XXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

PLOT HISTOGRAM OF LOOP OUTPUT ERRORS VS PHASE IN CEG. 

CONTRL=''1. 

IF(CONTRL.LT. 1Q.)G0TG 31 
CALL HIST1 {ERROR, SUM, SNR, VAR) 

31 CONTINUE 
CONTEL=0 . 

SNR=T0. +ALOG10 (A++2/(VABI*2.) ) 

XX (L) =SNR 
YY (L) -VAR 


Figure E-1 5, (Continucsd) 
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YY (L) =3UHT (VAR) 

5(LL)=;iNR 
STD (LL) =YY (L) 

LL=LL+1 

A=Af0.1 

29 FOanAT{lX,* INPUT SSE = • ,F7 ,3, 1 0X, • OUTPUT VARIANCE - SFQ.S) 

30 CONTINUE 
999 NPTSPP-100 

NPL0TS=1 
GO TO 4G 

CALL KMPLOT (XX,YY,NPTSPP,NPLOTS,XXL, YYL) 

A=Q.1 

40 CONTINUE 
GO TO 44 

44 CONTINUE 
SOGAIN=SOGAIN+. 125 

45 CONTINUE 

CALL NMPLOI (S , STP, 1 00, 6, XXL, YYL) 

STOP 
200 STOP 
END 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE QUANT (X , LEVELS) 

C ROUTINE QUANTIZES REAL X TO NUMBER OF LEVELS = LEVELS. 

G RETTURNS REAL X QUANTIZED, 

XX = FLOAT (LEVELS) 

X= (XX/J60.) *X 
I=INT(X) 

X=FLOAT(I) 

RETURN 

END 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE GH ANGE ( Y , INT VL, A ,SHA LLA , THE) 

INTEGER TME 

REAL Y(5000), INTVL 

TAU=INTVL 

E=EXr (-SMALLA’i'TAU) 

PSI= (A*A*S KALLA) /4 , 

CONS=SQKT(PSI*(1.- (R*S))) 

GALL ANORM(AN) 

Y(1) =AN«PSI 
DO 10 1=2, THE 
J=I-1 

CALL ANORM(AN) 

9 FORMAT (IX, 10E 10. 2) 

Y (I) =AN*CGNS+R*Y (J) 

1G CONTINUE 
RETURN 
END 

c XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE SMAPLL ( PHI ,L PGT 8 , PHI HEF, CONST, CTR2, SU MDIF) 

C ROUTINE PERFORMS DIGITAL FILTERING ON INCOMING OMEGA PHASE 

C MEASUREMENTS( PEL, BURST FILTER). OPERATION AND TIME CONSTANTS 

C CLOSELY RESEMBLE HARDWARE MAPLL. 

C *»CONST*' IS A VALUE TO BE SPECIFIED FOR MAX. LOOP CTF. VALUE 
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Figure E-15. (Cont’inued) 
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BEFORE ADJUSTING PHIfiEF. (LOOP TIME GOSSTAMT OF AnOUKT OF INTEGRA. 
B. SALTER, AVIONICS ENGINEERING CENTER, JULY 1976. 


C 


C 

G 


G 


G 

G 


C 


c 

c 

e 

c 

c 


10 


INTEGER PHI,PHIHEF,CONST,CrB2,LPCTR 
INTEGER DIFP,SOKDIF,HALFLN,LANE 
HALFLN=t)4 

COMPARf^INCOMING PHASE TO LOOP CONTROL WORD(LCW) {PHASE REFERENCE 


DIFF= PHI” PM JLREr 

KEEP A RUNNING SUN OF ABSOLUTE VALUES OF DIFFERENCES TO INDICATE 
SIGNAL VARIANCE. 

IF(DIFF).10, 50,20 

CHECK FOR EDGE OF LANE AMBIGUITY. 

IF(DIFF.LT.-HALFLN) GO TO 25 
SUMDIF=SUWDIF*IABS (DIFF) 


GO TO 40 

20 IF {DIFF, GT. HALFLN) GO TO 35 
SUKDIF=G0MDIFfIABS (DIFF) 

GO TO 30 

INCREMENT LOOP COUNTER. 

30 LPCTP=LPCT?.+ 1 

EXAMINE L)>CTH FOP IN CP. MORE THAN 
IF {LPCTP ,LE. CONST) RETURN 
LPCTR=0 

PilIHEF=PHIRBF+1 

CTR2=CTR2+1 

CHECK FOR MODULO 256 LANE COUNTS. 
IF(PHIHEF.GT.LANE) PHIREF=0 
RETURN 


•CONST* TIMES CORRELATED, 


DECREMENT LOOP COUNTER, 

40 LPCTP =LPCTR“ 1 

EXAMINE LPCTK FOP. DEGB. MORE THAN •CONST* TIMES COPRELATSD. 
IF (LPCTR ,GE, -CONST) RETURN' 

LPCTF=0 

PHIREF=PHIREF- 1 


GTE2=GTR2-1 

CHECK FOR MODULO 256 LANE COUNTS. 
IF ( P H IR E F . LT . 0) P HI R EF=L A N E 
50 RETURN 
25 Pai=P8I+LANS 

DIFF=PliI-PRIREF 
S 0 M D I F- S U M D I F + 1 A B S ( D 1 FF ) 

GO TO 30 
35 PHI=PHi-LANE 

DIFF=PHI~PHIBEF 
SUKDIF=SUMDI?+IABS (DIFF) 

GO TO 40 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXYXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

ROUTINE PLOTS DATA VECTOR ” X” WITH STABS, ONE POINT PER LIKE ON A 
1-128 SCALE. 

IF VALUES OF "X" GAN FALL BETWEEN 128 - 255 , SET LL=1 AND ROUTINE 


Figure E-IS. (Continued) 


- 112 - 






C KILL tiALVE EACH VALUE BEFORE PLOTTING. OTHERWISE, IF VALllFS OF 

C ”X*' AHE ALL LESS THAN 128, ROUTINE WILL PLOT STRAI3HT VALUES IF 

C LL=0 

C MUST DIMENSION X = VECTOR OF VALUES TO BE PLOTTED. 

C H.W. LILLEY, AVIONICS, JANUARY, 1976 

C MODIFIED BY h. SALTER, AVIONICS, APRIL, 1976. 

C 

INTEGER SK/* */,ST/' **/#LIN (128) 

INTEGER X(1000) 

DO 10 1=1,128 
10 LIN(I)=DK 

DO 20 1=1,200 
IF (LL.EQ.O) GO TO 12 
MM = X (I) /2+1 
GO TO 13 


12 continue 

MM=X (I) 

13 CONTINUE 
LIN (MM) =3T 
PRINT 11, LIN 

11 FORMAT (1X, 120A1) 
LIN (HH) =BK 
20 CONTINUE 
RETURN 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX . 

SUBROUTINE GH A NGB ( Y , I NT V L, A ,SM ALL A, TME) 

INTEGER THE 

PEAL I (5000) , INTVL 

TAU=INTVL 

R=EXP (-SKALLA«TAU) 

PS1= (A»A*SHALLA) /4, 

CONS = SQRT (PSI« (1 (B^R) ) )" 

CALL BNORM(AN) 
y (1) = AN*PSI 
DO 10 1=2, TME 
J=I-1 

CALL DNOKM(AN) 

Y (I) =AN*C0NS+R’>Y (J) 

10 CONTINUE 
RETURN 


XXXXXX'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE PLTPHS (LQUANT,IQUANT) -rt,P 

PLTPHS PLOTS THE QUANTIZED PHASE VALUE OF THE INPUT PHASE {♦) AND THE 

OUTPUT PHASE (♦} 

INTEGER BLANK, STAS, PT, CROSS 

DATA BLANK/’ ’/, ST AB/» * V, PT/' - V# CROSS/’ +’ / 

DIMENSION LINE(128) 

DO 10 K=1,128 
LINS(K) =BLANK 
LINE (LQUANT+1) =STAR 
C IF IT IS DESIRED THAT ONLY THE OUTPUT PHASE BE PLOTTED# SET 

C LINE(IQUANT) = BLANK 


10 


Figure E^15. (Conlinued) 
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OMGINAIi P^GE is 
OB) POOE quality 
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LINEdvfUANT+l) =CROSS 

IF<UJUANT.Ea.lgUAKT) LIMB(LQUAHT+1) =PT 
PRINT INLINE 
1 F0RBAT(1X,128A1) 

RETURN 

END 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE HN V AB (GA , AV E, VAR) 

REAL GA{5000) 

COMPUTE MEAN AND VARIANCE OF B.V.*S RETURNED FROM RANDOM NUMBER 
GENERATOR AND LOW PASS FILTER, 

saM=o. 

DO 5 1=1,5000 
SUW = SUM4-GA (I) 

5 CONTINUE 
AVE=SUM/5000. 

SUM2=0. 

DO 6 1=1,5000 
S0M2=SU212+ (GA (I) -AVE) **2 

6 CONTINUE 
VAE=SUK2/5000 
PRINT 7,AVE,VAS 

7 FORMAT(lX,2E20.y) 

RETURN 
END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROIIIINE HIST1 (NOISE, NUM, SNR, VAR2) 

INTEGER NOISE (5000) , NH (5000) 

REAL YYY (10 1) ,XXX (101) , XXXL (12) ,YVYL(12) ,XXL (12) ,yVL (12) 

DATA XXL/' S’ , •!' ,»G* , *N* , * A* , »L* , V » * N' ,'0' #*I' »*S', 'EV 
DATA YYL/»P» , • L* , • L’ , ’ » , » S » , • T » , ' D * , • • , * D' , » E ♦ , ' V ’ , * •/ 

DATA XXXL/'N • , '0* , ' 1 • * , < P ' , ' H * , ' A • , • S * , * E * , * V 

DATAyYyL/'h*,'I',»S','T»,'0»,'G*,»R',»AV,'M»,* *,* 'f* V 

COMPUTE AND PLOT PHOBAEILITY DENSITY OF THE INPUT NOISE BEFORE 
IT IS ADDED TO THE SIGNAL PHASE, QUANTIZED, AND SENT TO SMAPLL. 
ICEC=-50 
DO 1 fl I = 1 , 1 0 1 
YYY (I)=0 
DO 8 J=1,NUM 

NN (J) =I NT (FLOAT (NOISE (J) ) /3 .6) 

IF(NN (J) .EQ.ICEC) YYY (I) =YYY (I) f 1. 

307 F0RMAT<1X,I5, F6. 1) 

8 CONTINUE 
XXX (I)=FLOAT (ICEC) 

ICEC=ICEG+1 

18 CONTINUE 
DO 19 IK=1, 101 
YYY (IK) =YYY (IK) /FLOAT (NUM) 

19 CONTINUE 
GO TO 20 

GALL MMPLOT (XXX, YYY, 101,1, XXXL,YYYL) 

20 CONTINUE 
ISUK=0 

DO 115 1=1, NUM 


Figure E-15, (ConHniJsd) 
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lSUfS-^ISOH*NN(I) 

115 CONTINUE 

AVE=FLOAT (ISUM) /FLOAT (NUH) 

SUB=0 

DO 116 1=1, NUM 
ERH=FLOAT (NN{I) ) -AVE 
S0H=SUM+EEH**2 

116 CONTINUE 
VAE2^SUM/FLOAT (NUH) 

RETURN 

PRINT 1 17 ,AVE,VAR2,SNB 

117 FORMAT(1X,'MEAS=»,P10.5,1GX,*VARIANCE=*,F10.5,10X,'SNR= 
PBTUBN 

END 


A 




Figure E-15. (Cohtinued) 
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F, OSP Mtcrocomputer Routines . 

1 . Main Routine. The "main" routine in this version of the OSP is 


simply on initialization and idling routine. This program initializes all the memory 
words needed by the OSP data processing subroutines and then simply waits to 
be interrupted by the microcomputer interface module with new Omega data, as 
flow charted in Figure F-1 and listed In Figure F-2. This main routine would be 
replaced by a more meaningful area navigation routine in a full-blown ONS. 

The navigation routine could perform navigation processing functions in the time 


between Omega data Interrupts. 

2. Interrupt Service Routine. When the main program Is interrupted 
by the microcomputer interface module with new Omega data, the interrupt 
service routine is implemented as flow charted in Figure F-3 and listed in Figure 
F-4. The service routine first loads the new Omega data and updates a time-keeping 
register. It then cheeks a system status word to see If auto-sync is complete. If 
not, the data is passed to the auto-sync subroutine. If auto-sync is complete the 
time-keeping register is consulted to determine whether the current data sample 
is to be treated as Omega signal Information or a transmission gap sample. If 


the sample falls in a time-slot the SMAPLL2 subroutine is called and the data 


used to update the estimate of phose for that time-slot. If the data is noise-gap 


information the time-keeping register Is again consulted to determine if this Is the 
first measurement taken during the current gap. If so, the dota output subroutine 


is called and the rate-aiding for the second order SMAPLL2 is performed. Other- 
wise, the new data is simply forgotten and control is returned to the main routine 
to wait for the next Interrupt. 
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Figure F-1. Main (Initialization) Routine Flow Chart, 
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TIMCTE+1 

27 

106D 

DO 

02 



3NE 

FOUR 

28 

106F 

E6 

OD 



ING 

TiMCTE 

29 

1071 

A5 

OF 


FOOS 

LDA 

SYN'CST 

30 

107 3 

G9 

03 



CMP 

=03 

31 

1075 

DO 

ID 



3NE 

ELE V EN 

32 

1077 

AS 

QE 


FIVE 

LDA 

TIMCTR^I 

33 

1079 

10 

31 



3PL 

TWELEV 

34 

107B 

46 

11 


SIX 

LSS 

PHI 

35 

1 G7D 

20 

00 

13 


JSH 

SMAPLL 

36 

1080 

A9 

7F 


SEVEN 

LDA 

1111111 

37 

1082 

25 

OE 



AND 

TIKGTR+1 

38 

1084 

A 6 

42 



LDX 

TSLOT 

39 

1086 

S5 

03 



GKP 

LENGf E-1 , X 

40 

1088 

DO 

04 



3NE 

NINE 

41 

108A 

A9 

00 


EIGHT 

LDA- 

=00 

42 

108G 

85 

OE 



STA 

TIMGTF+1 

43 

108E 

68 



NINE 

PLA 


44 

T08F 

A8 




TAY 


45 

1090 

68 




PLA 


46 

1091 

AA 




TAX 


47 

1092 

68 




PLA 


48 

1093 

40 




RTI 


49 

1094 

20 

00 

11 

ELEVEN 

JSR 

SYNC 

50 

1097 

A5 

0E 



LDA 

TIMCTR+I 

51 

1099 

G9 

E8 



CMP 

= $E8 

52 

109B 

DO 

OC 



3NE 

B11 


ij V* 




F'^4, Microcomputer Program OSERV Listing. 
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53 109D A5 OD 

54 109F C9 Q3 

55 1QA1 DQ Of! 

56 10A3 A9 00 

57 10A5 35 OE 

58 10A7 85 OD 

59 10A9 4C 3E 10 

60 10 AC C9 01 

61 1GAE DO 3D 

62 10 BO A6 42 

63 10B2 ^5 20 

64 10B4 ,VG 0 8 

66 10B6 18 

66 ie-B7 75 29 

67 1GB9 95 29 

63 lOBB 4C CC 10 

69 lOBE 49 FP 

70 10 CO 18 

71 10C1 69 01 

72 10C3 85 20 

73 10C5 B5 29 

74 10C7 30 

75 10C8 F,5 2 0 

76 10GA 95 29 

77 10CC 20 00 14 

7 8 1 0 CF B A 

79 1 ODG 1 8 

80 10D1 65 42 

81 1003 AA 

82 1004 A9 OiJ 

83 1 C 06 95 30 

84 1008 35 31 

85 10DA 85 IF 

86 1O0C 85 20 

87 1O0F e6 42 

88 10E0 A9 09 

89 10B2 C5 42 

9 0 1GE4 00 A8 

91 10E6 A9 01 

92 10E8 85 42 

93 10EA 4C 8E 10 

94 ICED C9 14 

95 10SF 00 90 

96 TO FI A9 80 

97 10F3 85 OE 

98 IGI'5 4C 8E 10 

99 

END BASS 2; 


EDA TIKCTE 
CMP ~i03 
BKE B11 
LDA =00 
ST A TIMCTE+1 
STA TIHCTB 
B11 JMP NINE 
TWELEV CNP =^01 

3NE NINETN 
THIRIN LGX TSLGT 
LDA CTP.2 
3NI LA3EL1 
CLC 

ADC rw\TF-1»X 
STA BATE-1,X 
JKP FOUKTN 
LABELl EOF =$FF 
CLC 

ADC =01 
STA CTR2 
LDA KATE-1, X 
SEC 

SRC CTS2 
STA 3ATE-1,X 
FOURTN JSE OlffPOT 
FIFTN TXA 
CLC 

ADC TSLOT 
TAX 

LDA =00 
STA 3N-2,X 
STA SN-1,X 
STA CTPI 
S T* A CT B2 

SISTN INC TSLOT 
SEVTN LDA “09 

CMP TSLOT 
BNE NINE 
EIGtlTN LDA =01 

STA TSLOT 
JKP NINE 
NINETN CMP =GAP 
BNE NINE 
TKENTY LDA =X80 

STA TIKCTF.+ 1 
JKP NINE 
END 


0 EE BOBS 


3. Automatic Synehronizqtlon SubrouHne. The auto-sync procedure 
is explained in the main text/ and the microcomputer program is flow charted in 
Figure F-5. The microcomputer language listing appears in Figure F-6. 

4. SMAPLL2 Tracking Subroutine. The memory-aided second order 
phase locked loop (SMAPLL2) operation is also functionally outlined in the text. 
More detai ls on the performanee which can be expected from this tracking can^ 
figuration is given in Appendix E. The SMAPLL2 procedure as implemented in 
the OSP microcomputer is charted In Figure F-7/ and the program listing appears 
in Figure F-8, 

5. Etata_ Qujput Subroutine . The data outpiJt operation Is described 
in the text and is flow charted in Figure F-9. The microcomputer Output sub- 
routine is listed in Figure F-^-lO. 



Figure F-5. Automatic Synchronization Subroutine Flow Chart. 
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Figure F-5. {Continued) 














END PASS 1: 

0 EF 

ROKS 



1 

0000 




ORG 

00 

2 

0000 



SYMGST 

EQO 

16 

3 

0000 



TEMP 

EQU 

16 

4 

0000 



PHI 

EQtJ 

17 

5 

0000 



SPHI 

EQO 

72 

6 

0000 



RALFLK 

EQO 

63 

7 

0000 



LANE 

EQO 

127 

8 

0000 



FIRST 

EQO 

19 

9 

0000 



PREVPH 

EQO 

18 

10 

0000 



THRESH 

EQO 

05 

11 

0000 



KBYTE 

EQO 

125 

12 

0000 



BITTEN 

EQO 

90 

13 

0000 



TI!1CTR 

EQO 

13 

14 

0000 



PTR2 

EQO 

20 

15 

0000 



PTR1 

EQO 

69 

16 

0000 



PTR3 

EQO 

67 

17 

00 0 0 



PTR4 

EQO 

6 8 

18 

oooo 



LENGTH 

EQO 

04 

19 

0000 



STAST 

EQO 

03 

20 

0000 



STAST2 

EQO 

74 

21 

0000 



SHFCTR 

EQO 

21 

22 

0000 



NSE1IFT 

SQU 

23 

23 

0000 



CORCTR 

EQO 

25 

24 

0000 



MAXCOR 

EQO 

27 

25 

0000 



SYNCPT 

EQO 

29 

26 

1100 




ORG 

$1100 

27 

1100 

48 


SYNG 

PHA 


28 

1101 

98 



TYA 


29 

1102 

48 



PHA 


30 

1103 

8A 



TXA 


31 

1104 

48 



PHA 


32 

1105 

A 5 

OF 

TWO 

LDA 

SYNCST 

33 

1107 

DO 

5B 


3NE 

EIGHTH 

34 

1109 

A5 

1 1 

THREE 

LDA 

PHI 

35 

110B 

85 

48 


ST A 

SPHI 

36 

1 10D 

C9 

3F 


g:^p 

=HALFLN 

37 

110F 

90 

05 


3CC 

FIVE 

30 

1111 

A 9 

7f 

FOUR 

LDA 

= LANE 

39 

1113 

38 



SEC 


40 

1114 

E5 

48 


SBC 

SPHI 

41 

1 1 16 

A4 

OE 

FIVE 

LDY 

TinCTF«-1 

42 

1118 

CO 

01 


CPY 

=01 

43 

11 1A 

DO 

OA 


BNE 

EIGHT 

44 

11 1C 

85 

13 

SIX 

ST A 

FIRST 

45 

11 IE 

85 

12 


5TA 

PREVPH 

46 

11 20 

68 


SEVEN 

PL A 


47 

1121 

A8 



TAY 


48 

1122 

68 



PLA 


49 

1123 

A A 



TAX 


50 

1124 

68 



PLA 


51 

1125 

60 



RTS 


52 

1 1 26 

85 

48 

EIGHT 

ST A 

SPHI 


Figure F-6. Microcomputer Program SYNG Listing. 
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NINE 


53 

1128 

38 


54 

1129 

F.5 

12 

55 

112B 

BO 

05 

56 

1128 

49 

FF 

57 

112F 

18 


58 

1130 

69 

01 

59 

1132 

A4 

48 

60 

1134 

84 

12 

61 

1136 

EA 


62 

1137 

HA 


63 

1138 

EA 


64 

1139 

C9 

05 

65 

1 13B 

90 

04 

66 

113D 

18 


67 

11 3E 

4C 

42 11 

68 

1141 

38 


69 

1142 

A2 

7D 

70 

1144 

36 

59 

71 

1146 

CA 


72 

1147 

DO 

FE 

73 

1149 

A5 

OD 

74 

114B 

G9 

03 

75 

114D 

DO 

D1 

76 

1 14F 

A5 

OE 

77 

1151 

C9 

E8 

78 

1153 

DO 

CB 

79 

1155 

A5 

10 

80 

1157 

DO 

15 

81 

1159 

A9 

01 


82 115E 85 TO 

83 115B AS 13 

84 115F 85 48 

85 1161 4C 2a 11 

86 1164 G9 01 

87 1166 DO 03 

88 1168 4C 20 11 

89 116B 4C 06 12 

90 116E E6 OF 

91 1170 A9 00 

92 1172 85 OE 

93 1174 05 GD 

94 1176 A6 14 

95 1 178 135 04 

96 117a 85 15 

97 117C 26 4A 

98 117E 90 06 

99 1180 A9 01 

100 1182 35 10 

101 1 184 10 04 

102 1186 A9 00 

103 1188 85 10 

104 118A A6 43 

105 11 8G 05 5 A 

106 118E A4 44 

107 1190 FO 04 


FEC 

BBC PBFVPH 
8CS TEN 
EOR =$FF 
CLC 

ADC =01 
TEN LDY SPHI 

STY P3EVPH 
ELEVEN NOP 
NOP 
NOP 

THELEV CMP =THEESH 
3CC SIGNAL 
TBIPTN CLC 

JNP STORE 
SIGNAL SEC 
STORE LDX =NBYTE 
ROTATE aOL BITTEN-1»X 
DEX 

BNE aOIATF 
FOIJRTN LDA TIMCTF. 

CMP =$03 
BNE SEVEN 
LDA TIMCTF+1 
CMP =$E8 
BNE SEVEN 
FIFTN LDA TEMP 

BNE NINETN 
SIXTN LDA =01 
STA TEMP 

SEVTN LDA FIRST 
STA SPHI 
JMP NINE 
EIGHTN CMP =01 
BNE A19 
JMP SEVEN 
A19 JMP TWENTO 
NINETN INC SYNCST 
LDA =00 
STA TIMCTF+1 
STA TIMCTR 
TKENTA LDX PTR2 

LDA LENGTH, X 
TKENTB STA SHFCTR 
TKENTE 30L STAST2 
BCC PLDS 
LDA =01 
STA TEMP 
BPL TWENTF 
PLUS LDA =00 
STA TEMP 
TWENTF LDX PTES 

LDA BITTBN,X 
TKENTG LDY PTR4 
B20G BEQ T WE NTH 


F-6. (Continued) 
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1G8 

1192 

2A 


109 

1193 

88 


110 

1194 

10 

FA 

111 

1196 

2A 


112 

1197 

A 9 

00 

113 

1199 

65 

10 

114 

119b 

C9 

01 

115 

119D 

rO 

06 

116 

119F 

E6 

1A 

117 

11A1 

DO 

02 

118 

11A3 

E6 

19 

119 

11A5 

A4 

44 

12.0 

1 1A7 

C8 


121 

11A8 

CO 

08 

122 

1 1AA 

DO 

04 

123 

1 1 AC 

AO 

00 

124 

11AE 

E6 

43 

125 

11B0 

8U 

44 

126 

11B2 

C6 

15 

127 

11B4 

DO 

D4 

128 

11B6 

A5 

45 

129 

1 1B8 

DO 

08 

130 

1 1BA 

A9 

14 

131 

1 1BC 

85 

15 

132 

11BE 

Eb 

45 

133 

1 ICO 

10 

C4 

134 

1 1C2 

A5 

14 

135 

11G4 

C9 

07 

136 

11C6 

FO 

08 

137 

1 1C8 

E6 

m 

138 

11CA 

A9 

00 

139 

11CC 

85 

45 

140 

11CF 

FO 

A6 , 

141 

11 DO 

A9 

00 

142 

11D2 

85 

45 

143 

11D4 

85 

14 

144 

11D6 

85 

43 

145 

11D6 

85 

44 

146 

1 1DA 

A5 

19 

147 

11 DC 

C5 

IB 

148 

1 1 DE 

90 

18 

149 

11 EC 

DO 

06 

150 

11E2 

kb 

1 A 

161 

11E4 

C5 

1C 

152 

11E6 

90 

10 

153 

11E8 

A5 

19 

154 

1 1EA 

85 

1B 

155 

1 1EC 

A5 

U 

156 

11 EE 

85 

1C 

157 

11F0 

A5 

17 

158 

11F2 

85 

ID 

159 

1 1F4 

A 5 

18 

160 

11F6 

85 

IE 

161 

1 1F8 

A5 

18 

162 

UFA 

C9 

E8 


ROL A 
DEY 

BPL D20G 
TWENTH RGL A 
TKEWTI LDA =00 
ADC TEMT 
TWENTJ CMP =01 

BEQ TWEKTK 
INC CORCTR+I 
BNE TirfENTK 
INC CORCTE 
TWENTK LDY PTR4 
INY 

T«ENTL CPY =08 

BNE TWENTM 
LDY =00 
INC PTF3 
TWENTM STY PTR4 

DEC SHFCTF 
TKEKTN BNE TWFNTF 
TWENTO LDA PTPl 

BNE TWENTC- 
TKENTP LDA =20 

ST A S.'IFCTK 
INC PTR1 
BPL PLU^^ 
TWENTQ LDA' PTR2 
CMP =07 
BEQ T-JENTI 
TWENTT INC PTR2 
LDA =00 
STA PTR1 
BEQ T^IEKTA 
TWENT1 LDA =00 
STA PTR1 
STA PTR2 
STA PTR3 
STA PTB4 
LDA COBCTP 
CMP MAXCOR 
3CC TWENT4 
BNE T»<EKT2 
LDA CORCTE+1 
CMP MAXCQE+1 
BCG TMENT4 
TWENT2 LDA CORCTE 
STA MAXGOE 
LDA GORCTR+1 
STA MAXCOE+1 
TWENT3 LDA NSHIFT 
STA SYNCPT 
LDA NSHIFT+l 
STA SYNCPT* 1 
TWENT4 LDA KSHTFT+1 
CMP =$E8 


Figure F—6, (Goritinued) 
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163 

11FC DO ID 

164 

11FL A5 17 

165 

120C C9 03 

166 

1202 DO 17 

167 

1204 E6 OF 

168 

1206 A5 OD 

169 

1208 C5 ID 

170 

12GA DO OC 

171 

120C A5 OE 

172 

120E G5 IE 

173 

1210 DO 06 

174 

1212 L6 OF 

175 

1214 A9 80 

176 

1216 85 OE 

177 

1218 4C 20 11 

178 

121E E6 18 

179 

121D DO 02 

180 

12 IF E6 17 

181 

1221 A2 7D 

182 

1223 A5 5A 

103 

1225 30 03 

184 

1227 18 

185 

1228 90 01 

186 

122A 38 

187 

122L 36 59 

188 

122D GA 

189 

122E DO FB 

190 

1230 A9 00 

191 

1232 85 19 

192 

1234 85 1A 

193 

1236 A5 03 

194 

1238 85 4A 

195 

123A 4C.76 11 

196 



END PASS 2: 0 ERBOnS 


BNE TWSHT5 
LDA NS HI FT 
C1P =£03 
BNE TWENT5 
TKENT7 INC SYNCST 
XKENT8 LDA TTMCTR 
C«P SYNCPT 
BNE A25 
LDA TIMCTE+1 
CMP SYNCPT+1 
BNE A25 

TWENT9 INC SYNCST 
THIHTY LDA =£80 

STA TIHCTP+1 
A25 3M? seven 

TWENT5 INC NSHIFT+1 
BNE TWENT6 
INC NSHIFT 
TWEKT6 LDX =NBYTE 
LDA BITTBK 
BMI MINOS 
GLG 

3CC DOIT 
MINOS SEC 

DOIT aOL BITTBN-1,X 
DEX 

BNE DOIT 
LDA =00 
STA COBCTR 
STA COKGTF+1 
LDA STAST 
STA STAST2 
JMP TWENTA 
END 


Figure F-6. 
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END PftSS 1 


0 ERRORS 


1 OCOO 

2 0000 

3 OOOC 

4 0000 

5 0000 

6 0000 

7 0000 

8 0000 

9 0000 

10 oeoo 
1 1 0000 
12 0000 
13 0000 

14 0000 

15 0000 

16 0000 

17 13C0 

18 1300 48 

19 1301 4K 

20 1302 48 

21 1303 98 

22 1304 48 

23 1305 A5 OE 

24 1307 C9 81 

25 1309 DO 25 

26 130B A6 42 


27 

130D 

E5 29 

28 

130F 

30 OG 

29 

1311 

4 A 

30 

1312 

4 A 

31 

13YR 

85 4 6 

32 

1315 

B5 21 

3 3 

1317 

18 

34 

1318 

65 4 6 

35 

131A 

4C 2E 

36 

131 D 

4 9 FF 

37 

1 3 1 F 

18 

38 

1320 

69 01 

39 

1322 

4A 

40 

1323 

4A 

41 

1324 

S5 46 

42 

13 26 

B5 21 

43 

1328 

38 

44 

1329 

E5 46 

45 

132B 

20 9 0 

46 

132E 

95 21 

47 

1330 

A6 4 2 

48 

133? 

B5 21 

49 

1334 

38 

50 

1335 

E5 11 

51 

1337 

FO 2 A 

52 

13 39 

BO 41 


ORG 00 
TIHCTF; EQU 13 
CTHI EQU 31 
GTR2 SOU 32 
BATE SOU 42 
LCK EQU 34 
SN EQU 50 

PHI EQU 17 
TSL07 SOU 68 
PBEVPH EQ'U 18 
TE?1P EQU 16 
TFKP2 SQU 70 
HALFLN EQU 32 
CONST EQU 08 
f^0D64 EQU $1390 
SNR EQU S13G0 
ORG $1300 
SMAPLL PHA 
TXm 
PH A 
TYA 
PHA 

b1 LDA TlMGTR+1 

CMP =$81 
3NE El 

Cl LUX TSLOT 

LDA RATE-1, X 
BHI LAB ELI 
LSB A 
LSR A 
STA TE?1P2 
LDA LCW-1,X 
CLG 

ADG TEHP? 

13 JH? Dl 

LA-BELl EOR =$FF 



GLC 

ADC 

= 01 


LSR 

A 


LSR 

A 


STA 

TEMP 2 


LDA 

LCW-1 ,X 


■SEC 

SBC 

TEHP2 

D1 

JSR 

MOD 6 4 


STA 

LCW-1,X 

El 

LDX 

TSLOT 

LDA 

LCM-1,X 


SEC 

SBC 

PRI 

Fl 

3EQ 

N1 


3GS 

SI 


Figure F-8. MiGrocomputer Program SMAPLL2 Lisfing 


53 

133D 

A5 

1 1 


HI 

LDA 

PHI 

54 

133D 

38 




SFC 

LCW- 1 ,X 

55 

1 33F 

P5 

21 



SBC 

56 

1340 

20 

CO 

13 

Gl 

JSF 

SNR 

57 

1343 

C9 

20 



CMP 

=HALFLN 

58 

1345 

90 

22 



BCC 

01 

59 

1347 

C6 

IF 


11 

DEC 

CTP.1 

60 

1349 

A9 

FF 



LDA 

= 255 

61 

134B 

45 

1? 


J1 

EOR 

CTR1 

62 

134D 

18 




CLC 

= 01 

63 

134E 

89 

01 



ADC 

64 

13 50 

C9 

08 



CMP 

=CONST 

65 

1352 

DO 

OF 



BNE 

N1 

66 

1354 

A9 

00 



LDA 

=00 

67 

1356 

85 

IF 



STA 

CTR 1 

68 

1358 

C6 

20 


K1 

DEC 

CTR2 

69 

135^ 

D6 

21 


LI 

DEC 

LCM-1 ,X 

70 

135C 

B5 

21 


CKMOD 

LDA 

LCW- 1,X 

71 

13 5E 

20 

90 

13 

Ml 

JSR 

MOD64 

72 

1361 

95 

21 



STA 

LCW-1 ,X 

73 

1363 

68 



N1 

PL A 


74 

1364 

AB 




TAY 


75 

1365 

68 




PL A 


76 

1366 

AA 




TAX 


77 

1367 

68 




PL A 


78 

1368 

60 




RTS 

GTE 1 

79 

1369 

F.6 

IF 


01 

INC 

80 

136B 

A9 

08 


P1 

LDA 

=CONST 

81 

136D 

G5 

IF 



CMP 

CTE1 

82 

136F 

DO 

F2 



BNF 

N1 

83 

1371 

A9 

00 



LDA 

=00 

84 

1373 

85 

IF 



STA 

CTR1 

85 

1375 

E6 

20 


Q1 

INC 

CTR2 

86 

1377 

F6 

21 


B1 

INC 

LCW- 1,X 

87 

1379 

4C 

5C 

13 


JMP 

GKMOD 

88 

137C 

20 

CO 

1 3 

SI 

JSR 

SNR 

89 

137f 

C9 

20 


T1 

CMP 

=HALFI.N 

90 

1381 

BO 

E6 



BCS 

01 

9 1 

1383 

4C 

47 

13 


JSP 

11 

92 






END 


END 

PASS 2 


0 

EE ROSS 




Figure F-8, (Gohtinued) 


- 132 - 


END PASS 1: 0 ERSORS 


1 

00 00 



ORG 

00 

2 

0000 


TEMP2 

EQU 

70 

3 

1390 



ORG 

$1390 

4 

1390 85 

46 

MGD64 

ST A 

TEHP2 

5 

1392 8A 



TXA 


6 

1393 48 



PH A 


7 

1394 98 



TYA 


8 

1395 48 



PH A 


8 

1396 AS 

4 6 

TKO 

LDA 

TEMP 2 

10 

1398 30 

1 1 


3NI 

SEVEN 

11 

139 A 29 

4 0 

THREE 

AND 

=%O1OC0OO( 

12 

139C FO 

06 


QEQ 

SIX 

13 

139E AS 

46 

FO0H 

LDA 

TEMP 2 

14 

13AC 29 

3f 


AND 

=%30 1111 1' 

15 

13A2 85 

46 

FIVE 

STA 

TEMP 2 

16 

13A4 68 


SIX 

PLA 


17 

13A5 A8 



TAY 


18 

13A6 68 



PLA 


19 

13A7 AA 



TAX 


20 

13A8 AS 

46 


LDA 

TEMP2 

21 

13AA 60 



3TS 


22 

13AB 49 

FF 

SEVEN 

SOR 

= $FF 

23 

13AD 85 

46 


STA 

TEKP2 

24 

13AF E6 

46 


INC 

TEMP2 

25 

1381 A9 

3F 


LDA 

=63 

26 

1383 3 8 



SEC 


27 

1384 E5 

46 


SBC 

TEHP2 

28 

13E6 4G 

A2 

13 

JKP 

FIVE 

29 




END 


END 

PASS 2: 

0 

ERRORS 




Figure F-8, (Continued) 


END PASS 1: 


0 ERRORS 
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